in 


>l 


/O^i 


NAVAL  POSTGRADUATE  SCHOOL 

Monterey,  California 


THESIS 

A  STATISTICAL  ANALYSIS  OF  MONTHLY 
MONTEREY  PENINSULA  AND  THE  CARMEL 
CENTRAL  CALIFORNIA 

RAINFALL  FOR 
VALLEY  IN 

by 

David  Frederick  Davis 
March  1981 

Thesis  Advisor:                 P. A.  Jacobs 

Approved  for  public  release;  distribution  unlimited 


T199165 


SECURITY   CLASSIFICATION  OF   THIS  PAGE  (Wttmtt  Dmim  Entorod) 


REPORT  DOCUMENTATION  PAGE 


i.    report  numTTR 


2.  GOVT  ACCESSION  NO 


4.     TITLE  (and  Subtlilm) 

A  Statistical  Analysis  of  Monthly  Rainfall 
for  Monterey  Peninsula  and  the  Carmel 
Valley  in  Central  California 


7.     AUTHORS 

David   Frederick    Davis 


»   PERFORMING  ORGANIZATION  NAME  ANO  ADDRESS 

Naval  Postgraduate  School 
Monterey,  California  93940 


II.     CONTROLLING  OFFICE  NAME   AND  AOORESS 

Naval  Postgraduate  School 
Monterey,  California 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


J      RECIPIENT'S  CATALOG  NUMBER 


5.     TYPE  OF   REPORT  A   PERIOD  COVERED 

Master's  Thesis; 
March  1981 


S.  PERFORMING  ORG.  REPORT  NUMBER 


t.  CONTRACT  OR  GRANT  NUMBERS 


10.  PROGRAM  ELEMENT.  PROJECT  TASK 
AREA  *  WORK  UNIT  NUMBERS  ' 


12.     REPORT   DATE 

March    19  81 


IS.     NUMBER  OF  PAGES 

207 


14.    MONITORING  AGENCY  NAME  *   AOORESVIf  dlllotont  Irom  Controlling  Otlleu) 


IS.     SECURITY  CLASS,   lot  ihla  report; 

Unclassified 


IS*.     DECLASSIFICATION/  DOWNGRADING 
SCHEDULE 


l».     DISTRIBUTION   STATEMENT  (ol  :hii  Kaport) 

Approved  for  public  release;  distribution  unlimited 


17.     DISTRIBUTION  STATEMENT  (ol  tha  abatrmct  tntarad  In  Block  20,  II  dlllaront  from  Kapon) 


I*.     SUPPLEMENTARY  NOTES 


19.    KEY  WORDS  (Contlnuo  on  roworao  tida  II  nocoaaaty  m%4  Identity  by  block  mambor) 

Monthly  Rainfall  Models  Logistic  Analysis 

California  Rainfall  Logistic  Time  Series 

Seasonal  Rainfall  Forecasting 
ARMA  Models 


20.    ABSTRACT  (Cantimn  on  rororao  tlda  It  nocoaaavy  trxd  Identity  by  block  mmbot) 

This  thesis  presents  a  statistical  analysis  of  the  monthly 
rainfall  for  the  Monterey  Peninsula  and  the  Carmel  Valley  in 
Central  California.   The  analysis  begins  with  the  simple  first- 
order  autoregressive  Markov  model,  which  is  found  to  be  weak. 
Next,  2x2  contingency  tables  are  used  to  identify  predictors, 
one  of  which  is  found  to  be  January  rainfall.   Finally, 
logistic  analysis  is  used  to  quantify  the  predictive  ability 


do  ,; 


FORM 
AN  73 


1473  EDITION  OF   1  NOV  «•  IS  OBSOLETE 

S/N    0 102-014- 6«01   ! 


SECURITY  CLASSIFICATION  OF  THIS  PAOE  (Whon  Dot*  Kntarod) 


|u>WTy   CL*ttO'C*TlQW   Of    Twit  »»««f»  —  fW««  «■»«—  <■ 

if  January. 

This  paper  attempts  to  analyze  rainfall  time  series  in  the 
[statistical  sense.   No  attempt  is  made  to  provide  a  physical 
explanation  of  the  findings  from  the  point  of  view  of  a 
meteorologist . 


D  1  jSS^O     U73  2 

'N    0102-014-6601  iteu«t*v  cwami^catioh  o'  ▼*••  »»«t>»»~  o«*  »■»•»•*» 


Approved  for  public  release;  distribution  unlimited 


A  Statistical  Analysis  of  Monthly  Rainfall  for 
Monterey  Peninsula  and  the.  Carmel  Valley  in 
Central  California 


by 


David  Frederick  Davis 

Captain,  United  States  Army 

B.S.,  Colorado  School  of  Mines,  19  72 


Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degree  of 


MASTER  OF  SCIENCE  IN  OPERATIONS  RESEARCH 


ABSTRACT 

This  thesis  presents  a  statistical  analysis  of  the 
monthly  rainfall  for  the  Monterey  Peninsula  and  the  Carmel 
Valley  in  Central  California.   The  analysis  begins  with  the 
simple  first-order  autoregressive  Markov  model,  which  is 
found  to  be  weak.   Next,  2x2  contingency  tables  are  used 
to  identify  predictors,  one  of  which  is  found  to  be  January 
rainfall.   Finally,  logistic  analysis  is  used  to  quantify 
the  predictive  ability  of  January. 

This  paper  attempts  to  analyze  rainfall  time  series 
in  the  statistical  sense.   No  attempt  is  made  to  provide 
a  physical  explanation  of  the  findings  from  the  point  of 
view  of  a  meteorologist. 
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I.       INTRODUCTION 

A.  THE  PROBLEM 

The  Monterey  Peninsula  Water  Management  District,  in 
Central  California  coastal  area  has  as  one  of  its  responsi- 
bilities the  duty  to  recommend  and/or  impose  water  rationing 
on  its  constituents.   To  do  this  in  a  rational  way  requires 
the  District  to  have  some  formula  for  predicting  future  water 
availability.   Although  the  techniques  of  modern  meteorology 
are  becoming  more  sophisticated  and  exact  there  is  still  the 
inability  to  make  good  long-range  predictions.   This  thesis 
analyzes  three  series  of  Monterey  County  monthly  rainfall 
data  by  purely  statistical  methodology  in  order  to  identify 
possible  predictive  formulas. 

B .  NOTATION 

Rainfall  will  be  denoted  as   R.     which  will  represent 

t  ,m 

inches  of  rain  recorded  for  the   t   year  and  the   m    month, 
The  year  to  be  used  is  the  California  Water  Year  which  begins 
in  October  and  ends  the  following  September.   Thus   R,  ,   is 
the  monthly  rainfall  for  October  of  year  ,  ' 1"   and   Rg  «   is 
the  monthly  rainfall  of  May  of  year  '6'. 

An  overstruck  bar  as  in  R. .  will  indicate  the  arithmetic 
average  of  a  variable;  in  this  case  it  is  the  arithmetic 
average  of  all  years  and  months  of  rainfall.   R.     is  the 
average  of  rainfall  over  the  years  for  month  m;   R.. 
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represents  the  yearly  average  for  year  t. 


C.   METHODS  OF  ANALYSIS 

Three  methods  were  used  to  analyze  the  data.   The  first 
method  was  to  model  the  series  using  autoregressive  moving 
averages  as  described  in  Box  and  Jenkins  [Ref.  1].   The 
second  was  to  use  2x2  contingency  tables  to  identify  possible 
predictors.   The  third  was  logistic  regression  to  quantify 
the  findings  of  the  2x2  contingency  table  analysis.   These 
three  methods  will  be  described  in  further  sections  of  this 
paper. 

1.   ARMA(p,q)  Models 

A  widely  used  approach  to  time  series  modeling  pro- 
posed by  Box  and  Jenkins  is  the  ARMA(p,q)  model.   This  model 
is  actually  a  joining  of  two  types  of  model,  the  autoregres- 
sive and  the  moving  average. 

In  the  noation  of  Box  and  Jenkins: 
let   {Z,  ,  t=l,2,...,n}   be  a  time  series,  then  an  ARMA(p,q) 
process  may  be  written  as: 

t    Yl   t-1  Yp  t-p     t     1  t-1         q  t-q 

the   {a   ,  t=l,2,...,n}  are  assumed  to  be  random  shocks 

distributed  as  independent  and  identically  distributed  (iid) 

random  variables  with  mean  zero  and  variance  a*   and 

a 

Z  ~    Z,  -  Z..    The  further  assumption  of  normality  is  also 
usually  made. 

For  purposes  of  this  paper,  a  mapping  of   R. 
into   Z  ,  r=12(t-l)+m  was  made,  and  an  ARMA  analysis  was 
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then  conducted  on  this  index  transformed  series.   This 
analysis  is  described  in  section  III. 

2 .  2x2  Table  Analysis 

In  the  validation  of  section  IV  it  is  found  that 
the  ARMA  model  is  not  very  successful  in  describing  the 
data.   In  section  V  the  data  is  analyzed  by  means  of  2x2 
contingency  tables.   These  tables  are  good  tools  for  explo- 
ratory data  analysis  in  that  they  provide  a  visual  display 
of  the  data.   Statistical  procedures  based  on  the  null 
hypothesis  of  independence  can  be  used  to  quantify  the 
departure  from  independence.   The  theory  of  2x2  tables, 
and  contingency  tables  in  general  may  be  found  in  Fleiss 
[Ref.  3],  Dixon  and  Massey  [Ref.  5],  Brownlee  [Ref.  6], 
and  Mood,  Graybill,  and  Boes  [Ref.  7]. 

For  this  paper,  the  contingency  table  approach  is 
used  to  identify  a  month  or  group  of  months  of  a  year  whose 
rainfall  can  serve  as  a  predictor  for  the  rainfall  during 
the  remaining  months  of  the  year.   One  predictor  that  was 
suggested  is  the  rainfall  in  the  month  of  January. 

3 .  Logistic  Analysis 

Once  a  predictor  is  tentatively  identified  it  becomes 
necessary  to  quantify  the  degree,  direction  and  accuracy  of 
the  predictor. 

A  logistic  analysis  is  conducted  by  dividing  the  data 
for  a  year  into  two  sets,  the  predictor  or  control  set,  and 
predictand  or  complement  set.   For  this  analysis,  the  predictor 
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is  the  logged  anomaly  of  January  rainfall  for  the  year; 
that  is,  if  X   denotes  the  predictor  or  control  for  year 
t ,  then 

Xt  -  to(Rt,4»  "  k   Jj_  £n(Rt,4>  J-2 

(The  logarithm  is  used  to  better  symmetrize  the  model.)   The 
complement  is  the  raw  anomaly  of  the  total  rainfall  for  the 
immediately  subsequent  eleven  months;  that  is,  if  Y   denotes 
the  complement  for  year  t,  then; 

,      12  3 

:t  =   (  ^  Rt,m  +   ^  Rt+1, 

m=5   '     m=l     ' 


12  3 

^  -   (  !  R.  .  +   I      m) 


,   N-l   12  3 

~  ~TPT  .X   (  lK    Rt,m  +  I.    Rt+l,m) 
t=l   m=5        m=l 

Finally,  the  data  are  transformed  into  a  binary  representa- 
tion, relative  to  zero  as; 

(  0   if   Y1   <   0 
Y^  =  {  t  1.4 


t         I    1   if   Y't 


>   0 

In  section  VI  the  model  fit  is 

a+3  x 

P(Y=l|x=x)  = —5 1.5 

1  +  ea+6  x 

Where   x   is  as  before  and  P(Y=l|x=x)  is  interpreted  as: 
"the  conditional  probability  that  the  subsequent  eleven  month 
total  rainfall  will  be  above  its  mean,  given  that  the  logged 
anomaly  of  January  rainfall  was  'x'" . 
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II.   THE  DATA 

A.   GENERAL 

Three  data  sets  were  used  for  this  analyis.   The  location 
at  which  these  data  sets  were  gathered  is  shown  in  Figure  1. 
As  the  figure  indicates,  two  of  the  data  sets  are  on  the 
Monterey  Peninsula  proper,  while  the  third  set,  SC,  represents 
the  Carmel  River  Watershed  at  the  San  Clemente  Dam. 

Although  data  exists  in  all  cases  to  the  present,  all 
three  sets  were  truncated  at  September  of  19  74.   The  remain- 
ing data,  up  through  September  of  1980  was  reserved  for 
validation  of  the  models  and  methodology. 

The  data  coordinates  are: 

Data  set  RN:    36°  35*  42"  North  Latitude 
121°  54'  43"  West  Longitude 

Data  set  FL:    36°  35'  30"  North  Latitude 
121°  56'  30"  West  Longitude 

Data  set  SC:    36°  26'  12"  North  Latitude 
121°  42'  30"  West  Longitude 
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Figure  1.   Location  of  rainfall  data  sets 
and  the  years  available. 
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B.   DATA  SET  RN 

Data  set  RN  consists  of  monthly  rainfall  amounts  gathered 
by  Professor  R.J.  Renard,  Cooperative  Observer  for  the  National 
Weather  Service  Climatological  Station,  Monterey,  California. 
The  data  set  begins  in  June  1951  and  currently  terminates  in 
September  1980.   As  was  stated  above,  the  analysis  was  con- 
ducted only  on  that  data  between  and  including  October  1951 
and  September  19  74. 

1.   Raw  Data 

Appendix  A  contains  a  listing  of  data  set  RN .   Figure 
2  shows  the  raw  data  set.   Month  1  is  October  1951,  month  148 
is  January  1964,  and  up  to  month  288  which  is  September  1974. 
As  can  be  seen  the  data  are  strongly  seasonal.   This  is  enough 
to  indicate  that  the  series,  as  stated,  is  highly  non- 
stationary. 

The  data  presented  so  far  deals  with  only  monthly 
data.   Next  to  be  considered  is  the  series  of  yearly  total 
rainfalls.   The  results  are  shown  in  Figure  3  (Yearly  total 
rainfall) ,  4  (Correlogram  of  yearly  rainfall) ,  and  Table  1 
(Estimated  Autocorrelations).   In  this  case,  the  correlogram 
indicates  stationarity  and  independence  of  the  yearly  series. 
A  plot  of  the  lag  one  relationships,  Figure  5,  reinforces 
this  indication  of  independence. 
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Figure   2.      Monthly   rainfall   in   inches    for   data   set   RN 
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The  correlograms  and  Partial  Correlograms  to  follow 
indicate  the  95%  approximate  significance  levels  using 
dashed  lines.   For  development  of  these  significance 
levels  see  Box  and  Jenkins  [Ref.  1]. 
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Figure  3.   Yearly  total  rainfall  for  data  set  RN 
(1951  -  1974) . 
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Figure  4.   Correlogram  of  Yearly  total  rainfall 
for  data  set  RN. 


TABLE  1 

ESTIMATED  AUTOCORRELATIONS  OF  YEARLY 
TOTAL  RAINFALL  FOR  DATA  SET  RN 


LAG 

1 
2 
3 
4 
5 


AUTO  CO  RRE  L AT I ON  S 
VALUE        LAG 


VALUE 


200 

7 

.044 

109 

8 

-.109 

295 

9 

-.009 

186 

10 

-.139 

249 

11 

.248 

12 

-.211 
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Figure  5.   Lag  one  plot  of  yearly  rainfall  data 
for  data  set  RN . 


2 .   Swept  Data 

Pierce  [Ref.  9]  and  Hipel  [Ref.  11]  suggest  various 
ways  to  remove  the  seasonality  of  data  sets  like  RN,  FL,  and 
SC.   The  basic,  and  most  straight  forward  of  these  methods 
is  to  remove  the  various  monthly  means.   This  is  accomplished 
by  the  following  replacement: 
let 


R.»-  ™  —  R±-  ^  —R»  _ 
t  ,m    t ,m    m 

where   R.    represents  the  mean  of  the  month  m. 


II. 5 


One  statistic  that  is  a  byproduct  of  the  calculations 

of   R.    is   S?    defined  as  the  estimated  variance  of  the 

m        m 
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monthly  data  points: 


m   N-l    _,  t,m      m 

These  statistics  for  data  set  RN  are  shown  in  Table  2,  and 
illustrated  in  Figure  6 .   In  the  same  way  as  the  raw  data 
mapped  into  a  series,  a  series  is  created  from 


R.     as; 
t,m 


Zr  =  Rt     fr=12(t-l)+m  II. 7 


TABLE    2 
MONTHLY    MEANS    AND    VARIANCE    FOR    DATA    SET    RN 

MONTH  MEAN  VARIANCE 

1  .677  .5292 

2  2.478  3.6354 

3  3.355  5.9147 

4  4.124  5.4146 

5  2.592  4.5122 

6  2.639  3.7414 

7  1.708  2.7699 

8  .434  .2494 

9  .219  .1004 

10  .058  .0056 

11  .104  .0129 

12  .275  .3953 
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Figure  6.   Monthly  means  for  data  set  RN 
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Figure  7.   Monthly  rainfall  anomalies  in  inches 
for  data  set  RN 
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Figure  8.   Correlogram  of  the  monthly  rainfall 
anomalies  for  data  set  RN 


TABLE  3 

ESTIMATED  AUTOCORRELATIONS  OF  MONTHLY  RAINFALL 
ANOMALIES  FOR  DATA  SET  RN 
LAG      VALUE      LAG      VALUE 


1 

.249 

14 

-.057 

2 

-.090 

15 

-.109 

3 

-.059 

16 

-.006 

4 

.041 

17 

.067 

5 

.032 

18 

.004 

6 

-.035 

19 

.013 

7 

-.011 

20 

-.007 

8 

-.043 

21 

.063 

9 

-.073 

22 

.023 

10 

-.091 

23 

-.066 

11 

.076 

24 

-.044 

12 

-.020 

25 

.021 

13 

-.012 
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3 .   Logged  and  Swept  Data 

The  data  should  now  be  stationary  in  the  means. 
However,  as  seen  in  Table  2,  the  variances  of  monthly 
rainfall  amounts  are  not  homogeneous.   Kilmartin  [Ref.  10] 
discusses  various  transformations  of  the  data  to  remove 
this  heteroskedacity .   A  plot  of  the  variance  versus  mean, 
Figure  9  below,  indicates  that  the  logarithmic  transform 
of  the  data  might  be  useful. 
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Figure  9.   Plot  of  monthly  variance  against 
monthly  means  for  data  set  RN 
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Since  the  data  contain  zeros,  the  following  modified 
logarithmic  transformation  is  done 

N 


R 


■  m   =  ln(R.         +1)  -  i  I      £n(R.  m  +1)  II. 

t  ,m       t  ,m       N   ~,      t  ,m 


where  the  effect  of  the  addition  of  the  one  is  mostly  to 

preserve  the  mapping  of  zeros  into  zeros.   A  more  in  depth 

discussion  of  this  transformation  is  found  in  Kilmartin. 

The  mapping  is  performed  again  as  before  and   R'     and 

S1     are  calculated  in  a  manner  similar  to  II. 6  and 
.m 

shown  in  Table  4  and  Figures  10  and  11. 

TABLE  4 

MONTHLY  MEANS  AND  VARIANCE  FOR 
LOGGED  DATA  SET  RN 

MONTH 

1 

2 

3 

4 

5 

6 

7 

8 

9 
10 
11 
12 


MEAN 

VARIANCE 

.438 

.1549 

1.092 

.3454 

1.312 

.3538 

1.539 

.2003 

1.104 

.3795 

1.133 

.3706 

.854 

.2728 

.319 

.0767 

.176 

.0401 

.054 

.0045 

.094 

.0094 

.185 

.0849 
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Figure    10.      Monthly  means   of    logged   data   set   RN 
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Figure  11.   Plot  of  monthly  variance  against  monthly 
means  for  logged  data  set  RN 


37 


These  transformations,  the  logarithm  followed  by 
the  removal  of  the  monthly  means  of  the  logged  data,  result 
in  the  series  listed  in  Appendix  A  and  described  in  Figures 
12  and  13  with  Table  5. 

These  displays  indicate  that  a  suitably  stationary 
series  has  been  obtained.   Other  methods,  such  as  differencing, 
scaling,  and  Box-Cox  transformations,  see  Hipel  [Ref.  11], 
were  tried  but  with  less  success. 
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Figure  12a.   Logged  anomalies  of  monthly 
rainfall  for  data  set  RN.   Months  1  -148 


DRTfi    IS    FROM    NEHPN 


1  .233 


-1 . 103 


U& 


>(■>■::■"  i"";« 


II 


i!  :    J 


m " 


*.■>'■ .  i ""..■  i- 


•(  i 


i  / 


i  i       i  j 


;     |    I  '    Ml 


143 


MONTHS  OF    DRTR 


236 


Figure  12b.  Logged  anomalies  of  monthly 
rainfall  for  data  set  RN 
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Figure  13.   Correlogram  of  logged  anomalies  of 
monthly  rainfall  from  data  set  RN 
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TABLE  5 

ESTIMATED  AUTOCORRELATIONS  OF  LOGGED 
ANOMALIES  OF  MONTHLY  RAINFALL  FROM 
DATA  SET  RN 


AUTOCORRELATIONS 
LAG       VALUE        LAG        VALUE 


1 

.191 

14 

-.033 

2 

-.102 

15 

-.084 

3 

-.095 

16 

-.024 

4 

.071 

17 

.062 

5 

.056 

18 

-.015 

6 

-.053 

19 

-.008 

7 

-.009 

20 

.013 

8 

-.014 

21 

.013 

9 

-.022 

22 

-.038 

10 

-.069 

23 

-.024 

11 

.024 

24 

-.012 

12 

-.004 

25 

.032 

13 

.032 

C.   DATA  SET  FL 

The  label  for  these  data  derives  from  its  location, 
Forest  Lake,  on  the  Monterey  Peninsula,  in  Pebble  Beach, 
California.   Data  set  FL  consists  of  monthly  rainfall 
figures  gathered  by  the  California-American  Water  Company 
since  1896.   Although  this  data  set  started  quite  early, 
the  data  prior  to  19  37  has  frequent  missing  observations. 
Therefore,  this  data  set  is  taken  as  October  19  37  through 
September  1974,  with  October  1974  through  September  1980 
reserved  for  validation. 

Analysis  of  this  data  set  is  identical  to  that  of  data 
set  RN,  therefore  only  the  pertinent,  figures  and  tables 
are  shown. 
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Figure    14a.      Months    1    -    296    of    rainfall    in 
inches    for   data    set   FL 
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Figure  14b.   Months  297-444  of  rainfall 
in  inches  for  data  set  FL 
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Figure  15  .   Yearly  total  rainfall  for  data  set  FL 
(1937  -  1974) . 
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TABLE    6 
nNS    OF   ^ABLY   TOTAL   BMWM* 
ESTIMATBD  AUTOCORBB^XO^S    0^   ^ 

AUTOCORRELATIONS 

T 

LAG  VAME 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 


LAG 


.010 
.099 
.022 
.207 
.214 
.028 
-.003 
-.061 
-.138 
-.060 
.236 

-.119 
-.151 


14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 


VALUE 


-.028 

.105     . 
-.087 
.085 
-.018 
-.245 
-.237 
-.035 
-.031 
-.106 
-.096 
.120 
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2 .       Swept   Data 

TABLE    7 
MONTHLY    MEANS    AND   VARIANCE    FOR    DATA   SET    FL 


MONTH 

1 

2 

3 

4 

5 

6 

7 

8 

9 
10 
11 
12 


MEAN 

VARIANCE 

.744 

.4895 

2.235 

3.7444 

3.049 

4.2480 

3.537 

4.4240 

2.999 

5.9492 

2.724 

3.1743 

1.559 

2.4505 

.449 

.2033 

.153 

.0417 

.077 

.0081 

.115 

.0081 

.189 

.1350 
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Figure  17.   Monthly  means  for  data  set  FL 
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Figure  18a.   Months  1  -  296  of  rainfall  anomalies 
in  inches  for  data  set  FL 
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Figure  18b 


Months  297-444  of  rainfall  anomalies 
in  inches  for  data  set  FL 


CORRELOGRAM 


u 
a 

o 
u 
o 

3 
•X 


1  .0 


0.5 


0.0 


-0.5 


-1.0 


^  .  .-<~i 


25 


LAG 


Figure  19.   Correlogram  of  monthly  rainfall 
anomalies  for  data  set  FL 
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TABLE  8 

ESTIMATED  AUTOCORRELATIONS  OF  xMONTHLY 
RAINFALL  ANOMALIES  FOR  DATA  SET  FL 


AUTOCORRELATIONS 


LAG 


VALUE 


LAG 


VALUE 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 


244 

14 

-.053 

007 

15 

-.043 

056 

16 

-.002 

027 

17 

.039 

016 

18 

-.014 

026 

19 

-.011 

022 

20 

-.031 

021 

21 

.023 

051 

22 

.067 

020 

23 

-.007 

077 

24 

.041 

059 

25 

.039 

068 

3.   Logged  and  Swept  Data 
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Figure    20. 


Plot  of  monthly  variance  against  monthly 
means  for  data  set  FL 
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TABLE    9 
MONTHLY    MEANS    AND   VARIANCE    FOR   LOGGED    DATA    SET    FL 


MONTH 

1 

2 

3 

4 

5 

6 

7 

8 

9 
10 
11 
12 


MEAN 

VARIANCE 

.484 

.1440 

1.007 

.  3440 

1.269 

.2781 

1.398 

.2579 

1.218 

.3421 

1.184 

.3006 

.796 

.2717 

.338 

.0589 

.130 

.0234 

.071 

.0061 

.106 

.0064 

.146 

.0044 
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Figure  21.   Monthly  means  of  logged 
data  set  FL 
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Figure  22.   Plot  of  monthly  variance  against  monthly 
means  for  logged  data  set  FL 
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Figure  2  3a, 


Months  1  -  148  of  logged  rainfall 
anomalies  for  data  set  FL 
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Figure  2  3b. 


Months  149  -  444  of  logged  rainfall 
anomalies  from  data  set  FL 
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Figure  24. 


Correlogram  of  logged  anomalies  of  monthly 
rainfall  from  data  set  FL. 


TABLE  10 

ESTIMATED  AUTOCORRELATIONS  OF  LOGGED 
ANOMALIES  FROM  MONTHLY  RAINFALL  OF 
DATA  SET  FL 

AUTOCORRELATIONS 


LAG 

VALUE 

LAG 

VALUE 

1 

.185 

14 

-.052 

2 

-.020 

15 

-.021 

3 

-.081 

16 

-.020 

4 

.046 

17 

.040 

5 

.043 

18 

-.040 

6 

-.050 

19 

-.025 

7 

-.024 

20 

-.014 

8 

-.010 

21 

.007 

9 

-.027 

22 

.019 

10 

.004 

23 

-.015 

11 

.031 

24 

.076 

12 

.084 

25 

.047 

13 

.074 
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D.   DATA  SET  SC 

The  label  for  this  data  derives  for  its  location,  San 
Clemente  Dam,  on  the  Carmel  River  in  Central  Calif ronia, 
approximately  26  kilometers  southeast  of  data  sets  RN  and 
FL  on  the  Monterey  Peninsula.   Data  set  SC  consists  of 
monthly  rainfall  figures  gathered  by  the  California-American 
Water  Company  since  1926. 

Analysis  of  this  data  set  is  again  very  close  to  that 
of  the  previous  data  sets  and  only  the  displays  will  be 
given. 

1.   Raw  Data 


DRTfl  IS  FROM  NEMSC 


15.  160 


in 
UJ 

X 

u 


0.000  L 


I 


MONTHS    OF    DRTfl 


.48 


Figure  25a.   Months  1  (October  1926)  -  148 
(January  19  38)  of  rainfall  in  inches  for 
data  set  SC. 
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Fiugre  25b.   Months  149  -  444  of  rainfall 
for  data  set  SC 
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15.  ISO 
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Figure  25c.   Months  445  -  576  of  rainfall  for 
data  set  SC 
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Figure   26. 


Yearly  total  rainfall  for  data 
set  SC  (1926  -  1974) 
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Figure  27. 


Correlogram  of  yearly  total 
rainfall  for  data  set  SC 


TABLE  11 

ESTIMATED  AUTOCORRELATIONS  OF  YEARLY 
TOTAL  RAINFALL  FOR  DATA  SET  SC 

AUTOCORRELATIONS 


LAG 


VALUE 


LAG 


VALUE 


1 

-.050 

14 

-.109 

2 

.135 

15 

.161 

3 

-.158 

16 

.077 

4 

.168 

17 

.116 

5 

-.042 

18 

.025 

6 

.081 

19 

-.214 

7 

-.084 

20 

.066 

8 

-.217 

21 

.019 

9 

-.111 

22 

-.028 

10 

-.107 

23 

-.101 

11 

.158 

24 

-.030 

12 

-.260 

25 

-.042 

13 

-.110 
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2.   Swept  Data 


TABLE  12 

MONTHLY  MEANS  AND  VARIANCES 
FOR  DATA  SET  SC 


MONTH 

1 

2 

3 

4 

5 

6 

7 

8 

9 
10 
11 
12 


MEAN 

VARIANCE 

.698 

.5945 

2.175 

4.2382 

3. 940 

10.1783 

4.599 

8.4899 

4.353 

13.3443 

3.080 

5.2744 

1.700 

3.4486 

.431 

.1912 

.111 

.0451 

.017 

.0055 

.037 

.0125 

.103 

.1040 
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Figure  28.   Monthly  means  for  data  set  SC 
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Figure  29a.   Months  1-296  anomalies  in  inches 
for  data  set  SC 
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Figure  29b.   Months  29  7  -  5  76  anomalies  in 
inches  for  data  set  SC 
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Figure  30.   Correlogram  of  monthly  rainfall 
anomalies  for  data  set  SC 


TABLE  13 

ESTIMATED  AUTOCORRELATIONS  OF  MONTHLY 
RAINFALL  ANOMALIES  FOR  DATA  SET  SC 

AUTOCORRELATIONS 


LAG 


VALUE 


LAG 


VALUE 


1 

.140 

14 

-.088 

2 

-.039 

15 

-.051 

3 

-.021 

16 

-.006 

4 

.012 

17 

.015 

5   ' 

-.001 

18 

-.006 

6 

-.019 

19 

-.008 

7 

.003 

20 

-.027 

8 

-.013 

21 

-.003 

9 

-.038 

22 

-.122 

10 

-.038 

23 

-.006 

11 

.014 

24 

.011 

12 

.051 

25 

.011 

13 

.102 
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Logged  and  Swept  Data 
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Figure    31. 


Plot  of  monthly  variance  against  anomally 
means  for  data  set  SC. 


TABLE  14 

MONTHLY  MEANS  AND  VARIANCES  OF 
LOGGED  DATA  SET  SC 


ONTH 

MEAN 

VARIANCE 

1 

.444 

.1623 

2 

.961 

.3999 

3 

1.408 

.3949 

4 

1.583 

.3117 

5 

1.444 

.4928 

6 

1.243 

.3556 

7 

.817 

.3247 

8 

.320 

.0726 

9 

.092 

.0227 

10 

.015 

.0040 

11 

.031 

.0083 

12 

.075 

.0351 
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Figure    32.      Monthly  means   of    logged  data   set   SC 
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Figure    33. 


Plot  of  monthly  variance  against  monthly 
means  for  logged  data  set  SC 


61 


DfiTfl    IS    fROM    NE.JSC 


1  .-434 


in 

u 


cr 

q; 

a 
u 

o 
o 
o 

_l 


-1 .534 


MONTHS    OF    DRTR 


148 


1  .434 


u 


a; 

Cl 

a 
u 

'j 
o 


Lp, 


-1 .534 


DRTfl  IS  FROM  NEUSC 


i   III! 


I      \,1!'!, 


1?  II II  1     '     l't    I    III    (I  !il    till  «     I 


149 


MONTHS    OF    DfiTR 


29S 


Figure  34a.   Months  1  -  296  of  logged  rainfall 
anomalies  from  data  set  SC 
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Figure  34b.   Months  297  -  576  of  logged  rainfall 
anomalies  from  data  set  SC 
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Figure    35. 


Correlogram  of  logged  anomalies  of 
monthly  rainfall  from  data  set  SC 


TABLE  15 

ESTIMATED  AUTOCORRELATION  OF  LOGGED 
ANOMALIES  OF  MONTHLY  RAINFALL  FROM 
DATA  SET  SC 

AUTOCORRELATIONS 


LAG 


VALUE 


LAG 


VALUE 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 


096 
065 
066 
038 
012 
061 
023 
001 
056 
021 
007 
091 
091 


14 

-.057 

15 

-.022 

16 

-.023 

17 

.021 

18 

-.017 

19 

-.005 

20 

-.008 

21 

-.011 

22 

.092 

23 

-.019 

24 

.050 

25 

.002 
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III.   FIRST  ORDER  MARKOV  MODEL 

A.   THEORY 

As    first    shown   by   equation    1.1,    the    general   ARMA(p,q) 
model   is: 

lt  =   *1    Vl  +'--  +  ^p    2t.p   +at   -    81   at_L    -...-6q    at_q        III.l 

The  development  and  discussion  of  this  type  of  model  is 
contained  in  detail  in  Box  and  Jenkin  [Ref.  1]  and  Nelson 
[Ref.  8].   The  modeling  process  is  a  three  fold  procedure. 
The  parts  are: 

(1)  Identification 

(2)  Estimation 

(3)  Diagnosis. 

Identification  is  conducted  using  the  correlogram  and 
a  plot  of  the  partial-autocorrelations  (or  partial  correlo- 
gram) .   The  partial  autocorrelations  are  related  to  the 
autocorrelations,  see  Box  and  Jenkins  [Ref.  1] ,  Nelson 
[Ref.  8],  or  Richards  and  Woodall  [Ref.  12].   These  partial 
autocorrelations  are  used  to  determine  the  order  of  the 
moving  average  process  much  like  the  autocorrelations  may 
be  used  to  determine  the  order  of  the  auto-regressive 
process. 

Once  the  autocorrelations  and  partial  autocorrelations 
have  been  found,  the  degree  of  the  ARMA  may  be  estimated  by 
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techniques  described  in  Box  and  Jenkins,  Nelson  or  Richards 
and  Woodall.   Each  of  the  data  sets,  once  logged  and  swept, 
indicated  that  the  most  probable  model  was  an  ARMA(1,0)  or 
AR(1)  or  more  commonly  a  first-order  autoregressive  Markov 
model.   This  model  is  simply; 

K   =  p*t-l  +  at  "1.2 

where  the  p  is  the  autocorrelation  of  lag  one.  Thus,  this 
model  indicates  that  any  persistence  in  the  data  are  condi- 
tionally independent  of  the  past  given  the  lag  one  value. 

Subsections  B,  C,  and  D  below  show  this  model  as  applied 
to  the  three  data  sets  of  interest.   The  residuals  of  the 
model    Z.  -  PZ.  ,    are  examined.   The  residuals  appear  to 
be  independent,  however,  they  do  not  appear  to  be  normally 
distributed;  for  example,  there  is  a  high  peak  around  zero. 
One  possible  reason  for  this  discrepancy  may  be  the  dichotomy 
of  winter  and  summer  rain  as  indicated  in  Tables  2,  4,  7,  9, 
12,  and  14.   The  existence  of  months  with  zero  rainfall  during 
the  summer  suggests  that  one  should  consider  the  summer,  when 
rain  is  sparse,  completely  separate  from  the  winter  when  rain 
is  more  abundant.   Therefore,  also  shown  in  the  subsections 
below  is  the  autoregressive  model  applied  to  the  winter  months 
only.   This  is  accomplished  by  stripping  out  months  9  through 
12  (June  through  September)  of  the  data  sets  and  treating  the 
remaining  data  as  a  continuous  set.   In  other  words,  the 
first  ten  months  are  then 

Rl,l'  Rl,2'  Rl,3/Rl,4/Rl,5'  R1,6'R1,8'R2,1'R2,2* 
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The  appropriate  correlograms  and  partial  correlograms  are 
displayed  prior  to  the  model  applications. 

B.   DATA  SET  RN 

1.   Twelve  Month  Series 

This  data  set  is  described  in  section  II. b.   The 
remaining  diagnostic  device  needed  is  the  partial  correlogram 
of  Figure  36  and  the  corresponding  values  in  Table  16. 
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Figure  36.   Partial  correlogram  of  the  logged 
rainfall  anomalies  of  data  set  RN 
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TABLE  16 

ESTIMATED  PARTIAL-AUTOCORRELATIONS  FOR 
LOGGED  RAINFALL  ANOMALIES  OF  DATA  SET  RN 

LAG        VALUE        LAG       VALUE 


1 

.191 

14 

-.040 

2 

-.144 

15 

-.073 

3 

-.047 

16 

.001 

4 

.092 

17 

.054 

5 

.006 

18 

-.067 

6 

-.058 

19 

.039 

7 

.037 

20 

.012 

8 

-.034 

21 

-.001 

9 

-.028 

22 

-.043 

10 

-.056 

23 

.013 

11 

.048 

24 

-.043 

12 

-.041 

25 

-.034 

13 

.047 

The  model  of  interest  is  then 

i't-  -isiit-i  +  at  XI1-3 

where  the  random  shocks   {a. }   are  assumed  to  be  distributed 

iid  N(0,a2)   and  a2    is  estimated  as 
a         a 

n^T  J  (zt  -  ■191zt-l)  IIX-4 


The  goodness  of  this  fit  may  be  viewed  in  two  ways. 
Firstly,  are  the  residuals,   {a. }   independent?   Secondly, 
are  the  residuals  distributed  as  Normal  (Gaussian)  random 
variables?   A  plot  of  the  residuals  follows  in  Figure  37 . 
The  question  of  independence  is  addressed  in  Figure  38 
(Correlogram) ,  Figure  39  (Lag  one  plot) ,  Figure  40  (Residuals 
vs.  lag  one) ,  and  Table  17  (Turning  points) .   For  a  discussion 
of  the  usefulness  of  the  turning  points  see  Kendall  [Ref.  14]. 
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All  of  these  displays  and  tests  tend  to  indicate  that  the 
residuals  are  in  fact  serially  independent.   The  statistics 
of  the  residuals  are  in  Table  18.   A  Normal  Plot  of  the 
residuals  (Figure  41) ,  in  which  the  sample  is  normalized  by 
removing  the  mean  and  scaling  by  the  standard  deviation  and 
then  plotted  on  normal  paper,  should  yield  a  nearly  straight 
line  corresponding  to  the  dashed  line  of  the  figure.   The 
Normal  Plot  accompanied  by  the  sample  histogram  (Figure  42) 
addresses  the  normality  of  these  data.   As  may  be  seen  from 
the  kurtosis,  the  fluctuations  of  the  sample  CDF  near  the 
midpoint,  and  the  peak  of  the  histogram,  the  normality  of 
this  data  are  questionable.   To  confirm  this  a  chi-squared 
goodness  of  fit  test  was  conducted  yielding  a  value  of  49.18 
with  17  degrees  of  freedom,  again  rejecting  any  hypothesis 
of  normality  at  a  significance  level  of  5x10 
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Figure  37.   First  order  Markov  residuals  from  logged 
rainfall  anomalies  of  data  set  RN .   Months  149  -  292 
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DATA  IS  FROM  RNLGSH 
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Figure  37b.   First  order  Markov  residuals  from 
logged  rainfall  anomalies  of  data  set  RN. 
Months  149  -  292 
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Figure  38.   Auto  correlations  of  residuals  from 
first  order  Markov  process  applied  to  the  logged 
rainfall  anomalies  of  data  set  RN 
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Figure  39.   Lag  one  plot  of  first  order  Markov 
residuals  from  logged  rainfall  anomalies  of 
data  set  RN 
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Figure  40.   First  order  Markov  residuals  versus 
lag  one  data  point  from  logged  rainfall  anomalies 
of  data  set  RN 


TABLE  17 

ACTUAL   AND   EXPECTED   NUMBER   OF    TURNING   POINTS 
AND   ACTUAL    AND   EXPECTED    PHASE    FREQUENCIES 
FOR   THE    FIRST    ORDER   MARKOV    RESIDUALS    FROM 

THE    LOGGED    RAINFALL   ANOMALIES    OF    DATA   SET    RN 

NUMBER   OF    TURNING   POINTS    =191 
E[P]    =    190.667  V[P]      =    15.899 

PHASE    LENGTHS 


1 
2 
3 
4 
5 
6 
7 


OBS. 

E[*] 

117 

118.8 

56 

52.1 

15 

14.9 

3 

3.2 

0 

.6 

0 

.1 

0 

0.0 
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TABLE  18 

GENERAL  STATISTICS  OF  FIRST  ORDER 
MARKOV  RESIDUALS  FROM  LOGGED  RAINFALL 
ANOMALIES  OF  DATA  SET  RN 


Moments 


Mean 
Variance 
Skewness 
Kurtosis 


-.001 
.117 

-.066 
.523 


Percentiles 


Minimum 

Lower  Sixteenth 

Lower  Eight 

Lower  Qua r tile 

Median 

Upper  Quartile 

Upper  Eight 

Upper  Sixteenth 

Maximum 


■1.141 

■  .745 

•  .463 

•  .174 

■  .014 
.211 
.436 
.706 

1.246 
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Figure  41.   Standardized  normal  plot  of  first 
order  Markov  residuals  from  logged  rainfall 
anamlies  of  data  set  RN 
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Figure  42.   Histogram  of  first  order  Markov  residuals 
from  logged  rainfall  anomalies  of  data  set  RN 


2 .   Winter  Series 

As  stated  above,  the  number  of  summer  months  with 
zeros  indicated  that  a  look  at  the  winter  months  only  might 
be  worthwhile.   The  Figures  4  3  (Winter  months) ,  44  (Correlo- 
gram) ,  and  45  (Partial  correlogram) ,  which  deal  only  with 
winter  months,  still  indicate  a  first  order  autoregressive 
Markov  model  as; 

Zt  =  .218Zt_1  +  at_x  III. 5 
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Figure  43.   Winter  months  only  of  logged  rainfall 
anomalies  of  data  set  RN 
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Figure  44.   Correlogram  of  winter  months  only  of 
logged  rainfall  anomalies  from  data  set  RN 
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Figure  45.   Partial  autocorrelations  of  winter 
months  only  logged  rainfall  anomalies  from  data 
set  RN 
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Now,  as  with  the  full  twelve  month  model,  a  look 
at  the  residuals  yields  the  Figures  46  (Residuals) ,  47 
(Correlogram) ,  and  48  (Lag  one  plot) ,  49  (Residuals  vs. 
lag  one) ,  and  Table  19  (Turning  points) .   It  appears  that 
the  residuals  are,  in  fact,  independent.   This  is  similar 
to  the  twelve  month  model. 

The  question  of  the  normality  of  the  residuals  is 
addressed  by  Table  19  and  Figures  50  (Normal  plot)  and  51 
(Histogram) .   The  results  of  these  plots  and  a  basic 
chi-squared  goodness  of  fit  of  22.11  with  17  degrees  of 
freedom  indicate  that  this  winter  month  data  set  is  much 
more  normal  than  was  its  twelve  month  counterpart.   This 
chi-squared  value  is  significant  at  the  .181  level. 
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Figure  46a.   First  order  Markov  residuals  of  logged 
rainfall  anomalies  for  winter  months  only  of  data 
set  RN.   Years  1-12 
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Figure  46b.  First  order  Markov  residuals  of  logged 
rainfall  anomalies  for  winter  months  only  of  data 
set  RN  .   Years  12  -  24 
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Figure  47.   Correlogram  of  first  order  Markov  residuals 
of  logged  rainfall  anomalies  for  winter  months  only  of 
data  set  RN 
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Figure  48.  Lag  one  plot  of  first  order  Markov  residuals 
from  logged  rainfall  anomalies  for  winter  months  only  of 
data  set  RN 
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Figure  49 .   First  order  Markov  residuals  versus  lag 
one  data  point  from  lagged  rainfall  anomalies  of 
winter  month  only  data  from  data  set  RN 


TABLE  19 

ACTUAL   AND    EXPECTED   NUMBER   OF    TURNING 
POINTS    AND   ACTUAL   AND    EXPECTED    PHASE 
FREQUENCIES    FROM   THE    FIRST   ORDER  MARKOV 
RESIDUALS    OF    THE    LOGGED    RAINFALL   ANOMALIES 
OF    DATA   SET    RN 

NUMBER   OF    TURNING   POINTS    =    129 
E[P]    =    126.667  V[P]       =    15.84 

PHASE    LENGTHS 


1 
2 
3 
4 
5 
6 
7 


OBS. 

E[*] 

82 

78.8 

38 

34.5 

7 

9.9 

1 

2.1 

1 

.4 

0 

0.0 

0 

0.0 
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TABLE  2  0 

GENERAL  STATISTICS  OF  FIRST  ORDER 
MARKOV  RESIDUALS  FROM  LOGGED  RAINFALL 
ANOMALIES  OF  DATA  SET  RN 


t 

Moments 

Mean 

-.001 

Variance 

.247 

Skewness 

-.161 

Kurtosis 

-.521 

Percentiles 


Maximum 

Lower  Sixteenth 

Lower  Eight 

Lower  Qurtile 

Median 

Upper  Quartile 

Upper  Eight 

Upper  Sixteenth 

Maximum 


■1.148 

•  .845 

•  .616 
■    .368 

.025 
.339 
.591 
.762 
1.078 
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Figure  50.   Standardized  normal  plot  of  first  order 
Markov  residuals  from  logged  rainfall  anomalies  of 
winter  months  only  from  data  set  RN 
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Figure  51.   Histogram  of  first  order  Markov  residuals 
from  logged  rainfall  anomalies  of  winter  month  only 
of  data  set  RN 


C.   DATA  SET  FL 

As  in  the  previous  section  on  the  data  sets,  the  analysis 
of  section  III.B  above  carries  forward  fairly  well  to  data 
sets  FL  and  SC.   This  section,  and  the  following,  contain 
only  the  Figures  and  Tables  corresponding  to  those  in  the 
previous  section  on  data  set  RN. 
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1.       Twelve   Month   Series 
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Figure    52. 


Partial  correlogram  of  the  logged 
rainfall  anomalies  for  data  set  FL 


TABLE  21 

ESTIMATED  PARTIAL  AUTOCORRELATIONS  FOR 
LOGGED  RAINFALL  ANOMALIES  OF  DATA  SET  FL 


LAG 

VALUE 

LAG 

VALUE 

1 

.185 

14 

-.072 

2 

-.057 

15 

.011 

3 

-.069 

16 

-.019 

4 

.077 

17 

.027 

5 

.016 

18 

-.048 

6 

-.068 

19 

.006 

7 

.010 

20 

-.007 

8 

-.008 

21 

.003 

9 

-.040 

22 

.018 

10 

.022 

23 

-.025 

11 

.031 

24 

.081 

12 

.067 

25 

.012 

13 

.055 

~  i 

Zt  " 

.185^^  +  at 

III. 6 
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Figure  5  3a.   First  order  Markov  residuals  from 
logged  rainfall  anomalies  of  data  set  FL. 
Months  1-296 
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Figure  5  3b.   First  order  Markov  residuals  from 
logged  rainfall  anomalies  of  data  set  FL. 
Months  29  7  -  444. 
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Figure  54.   Autocorrelations  of  residuals  from  first 
order  Markov  process  applied  to  the  logged  rainfall 
anomalies  of  data  set  FL 
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Figure  55.   Lag  one  plot  of  first  order  Markov 
residuals  from  logged  rainfall  anomalies  of 
data  set  FL 
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Figure  56.   First  order  Markov  residuals  versus 
lag  one  data  points  from  logged  rainfall  anomalies 
of  data  set  FL 
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TABLE  22 

ACTUAL  AND  EXPECTED  NUMBER  OF  TURNING 
POINTS  AND  ACTUAL  AND  EXPECTED  PHASE 
FREQUENCIES  FROM  THE  FIRST  ORDER  MARKOV 
RESIDUALS  FROM  DATA  SET  FL 

NUMBER  OF  TURNING  POINTS  =  294 
E[P]  =  294.667       V  [P  ]  =  15.93 

PHASE  LENGTHS 


D 

OBS. 

E[*] 

1 

188 

183.8 

2 

79 

80.7 

3 

19 

23.2 

4 

6 

5.0 

5 

0 

.9 

6 

1 

.1 

7 

1 

0.0 

8 

0 

0.0 

9 

0 

0.0 

10 

0 

0.0 

TOTALS 

294 
TABLE  2  3 

293.7 

GENERAL  STATISTICS  OF  FIRST  ORDER  MARKOV 

RESIDUALS  FROM  LOGGED  RAINFALL  ANOMALIES 

OF  DATA  SET  FL 

Moments 


Mean 

Variance 
Skewness 
Kurtosis 


.000 

.163 

-.044 

.717 


Percentiles 


Minimum 

Lower  Sixteenth 

Lower  Eight 

Lower  Quartile 

Median 

Upper  Quartile 

Upper  Eight 

Upper  Sixteenth 

Maximum 


•1.276 

■  .702 

■  .426 
•    .156 

■  .012 
.184 
.481 
.648 

1.124 
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Figure  57.   Standardized  normal  plot  of  first  order 
Markov  residuals  from  logged  rainfall  anomalies  of 
data  set  FL  
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Figure  58.   Histogram  of  first  order  Markov  residuals 
from  logged  rainfall  anomalies  of  data  set  FL 
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This  data  set  yielded  a  chi-square  value  of  1Q7.66 
for  17  degrees  of  freedom.   The  significance  of  the  value 
is  in  the  zero  plus  range. 
2.   Winter  Series 
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Figure  59a.   Years  1  -  12  of  winter  months  only  of 
logged  rainfall  anomalies  of  data  set  FL 
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Figure  59b.   Years  13  -  37  of  winter  months  only, 
logged  rainfall  anomalies  of  data  set  FL 
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Figure  60.   Correlogram  of  winter  months  only, 
logged  rainfall  anomalies  from  data  set  FL 
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Figure  61.   Partial  correlogram  of  winter  months  only, 
logged  rainfall  anomalies  from  data  set  FL 


These  displays  indicate  a  model  like 
z"  =  .199Z.  ,  +  a. . 
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Figure  62a.   Years  1-24,  first  order  Markov 
residuals  of  logged  rainfall  anomalies  for  winter 
months  only,  data  set  FL 
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DATA  IS  FROM  FLLGSW 


1.  13£ 


D 
Q 
t— t 

lo 

u 

> 

o 

a: 

u 
a 
ct 
o 


01 

a. 


-1  .274 


193 


WINTER  MONTHS  OF  DflTfl 


288 


Figure  62b.   Years  25  -  37,  first  order  Markov 
residuals  of  logged  rainfall  anomalies  for 
winter  months  only,  data  set  FL 
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Figure  63.   Correlogram  of  first  order  Markov 
residuals  of  lagged  rainfall  anomalies  from  winter 
months  only,  data  set  FL 
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Figure  64.   Lag  one  plot  of  first  order  Markov 
residuals  from  logged  rainfall  anomalies  of  winter 
months  only,  data  set  FL 
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Figure  65.   First  order  Markov  residuals  versus  lag 
one  data  point  from  logged  rainfall  anomalies  of 
winter  months  only,  data  set  FL 
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TABLE  2  4 

ACTUAL  AND  EXPECTED  NUMBER  OF  TURNING 
POINTS  AND  ACTUAL  AND  EXPECTED  PHASE 
FREQUENCIES  FROM  THE  FIRST  ORDER  MARKOV 
RESIDUALS  OF  THE  LOGGED  RAINFALL  ANOMALIES 
OF  DATA  SET  FL 

NUMBER  OF  TURNING  POINTS  =  209 
E[P]  =  196  V[P]  =  15.902 

PHASE  LENGTHS 


1 

OBS. 

E[*] 

1 

138 

122.1 

2 

59 

53.5 

3 

8 

15.4 

4 

2 

3.3 

5 

1 

.6 

6 

0 

.1 

7 

0 

0.0 

8 

0 

0.0 

TABLE  2  5 

GENERAL  STATISTICS  OF  FIRST  ORDER 

MARKOV  RESIDUALS  FROM  LOGGED  RAINFALL 

ANOMALIES  OF  WINTER  MONTHS  ONLY, 

DATA  SET  FL 

Moments 


Mean 

Variance 
Skewness 
Kurtosis 


.000 

.234 

-.079 

-.376 


Percentiles 


Minimum 

Lower  Sixteenth 

Lower  Eight 

Lower  Quartile 

Median 

Upper  Quartile 

Upper  Eight 

Upper  Sixteenth 

Maximum 


■1.274 

■  .798 

■  .570 

■  .315 
.011 
.323 
.551 
.748 

1.136 
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Figure  66.   Standardized  normal  plot  of  first 
order  Markov  residuals  from  logged  rainfall 
anomalies  of  winter  months  only,  data  set  FL 
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Figure  67.  Histogram  of  first  order  Markov 
residuals  from  logged  rainfall  anomalies  of 
winter  months  only,  data  set  FL 
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The  chi-squared  was  calculated  at  15.35  for  17 
degrees  of  freedom.   This  is  a  significance  level  of  .570, 
thus  indicating  possible  normality. 

D.   DATA  SET  SC 

1.   Twelve  Month  Series 
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Figure  68.   Partial  correlogram  of  the  logged 
anomalies  for  data  set  SC 
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TABLE  26 

ESTIMATED  PARTIAL  AUTOCORRELATIONS 
FOR  LOGGED  RAINFALL  ANOMALIES  OF 
DATA  SET  SC 


LAG 


VALUE 


LAG 


VALUE 


1 

.096 

14 

-.066 

2 

-.075 

15 

.004 

3 

-.053 

16 

-.022 

4 

.046 

17 

.007 

5 

-.004 

18 

-.014 

6 

-.061 

19 

.004 

7 

.042 

20 

-.019 

8 

-.016 

21 

.002 

9 

-.061 

22 

.106 

10 

-.000 

23 

-.046 

11 

-.001 

24 

.059 

12 

.079 

25 

-.015 

13 

.084 

This  information  yields  the  model  as 

Z,  =  .096Z.  ,  +  a, . 
t        t-1     t 
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Figure  69a.   First  order  Markov  residuals  from  logged 
rainfall  anomalies  of  data  set  SC.   Months  1  -  148. 
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Figure  6  9b.   First  order  Markov  residual  from 
logged  rainfall  anomalies  of  data  set  SC. 
Months  149  -  444 
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Figure  6  9c.   First  order  Markov  residual  from  logged 
rainfall  anomalies  of  data  set  SC.   Months  445  -  596, 
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Figure  70.   Autocorrelations  of  residuals  from  first 
order  Markov  process  applied  to  the  logged  rainfall 
anomalies  of  data  set  SC 
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Figure  71.   Lag  one  plot  of  first  order  Markov 
residuals  from  logged  rainfall  anomalies  of 
data  set  SC 


1  .344 


x 

X    x 
) 

;   x  » 

x    x 


c*x 


x,     x>x 

*       X* 


a 

(/) 

u 

a. 


-*«_ 


XX 

<* 

X 


-  X      x?*    A\ 


*x 
x    *■ 


„      xxx  x  x 
X   #X 


X*, 


*?£*    ■• 


*xx 

X     x 


V 


x  *     .. 
X 
X  .. 


*X 

#x 

i:  x 


-1.515 

-1.515 


X 

xXx 

.** 
x 


X  Sx 

XX  x 

X 


X      * 


1.344 


LAGGED    ELEMENTS 
H(t-l) 


Figure.  72.   First  order  Markov  residuals  versus 
lag  one  data  points  from  logged  rainfall  anomalies 
of  data  set  SC 
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TABLE  2  7 

ACTUAL  AND  EXPECTED  NUMBER  OF  TURNING 
POINTS  AND  ACTUAL  AND  EXPECTED  PHASE 
FREQUENCIES  FROM  THE  FIRST  ORDER  MARKOV 
RESIDUALS  OF  DATA  SET  SC 


NUMBER  OF  TURNING  POINTS  =36  7 
E[P]  =  382.667       V  [P  ]  =  15.9 


PHASE  LENGTHS 

D  OBS 

1  226 

2  95 

3  34 

4  8 

5  2 

6  1 

7  1 

8  0 

9  0 
10  0 

TOTALS  36  7 


E[*] 

238.8 

104.9 
30.1 
6.6 
1.2 
.2 
0.0 
0.0 
0.0 
0.0 

381.7 


TABLE  2  8 

GENERAL  STATISTICS  OF  FIRST  ORDER 

MARKOV  RESIDUALS  FROM  LOGGED 
RAINFALL  ANOMALIES  OF  DATA  SET  SC 

Moments 


Mean 

Variance 
Skewness 
Kurtosis 


000 
209 
057 
805 


Percentiles 


Minimum 

Lower  Sixteenth 

Lower  Eight 

Lower  Quartile 

Median 

Upper  Quartile 

Upper  Eight 

Upper  Sixteenth 

Maximum 


•1.515 

•  .744 

■  .462 

•  .205 

■  .030 
.207 
.506 
.783 

1.344 
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Figure  73.   Standardized  normal  plot  of  first  order 
Markov  residuals  from  logged  rainfall  anomalies  of 
data  set  SC 
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Figure  74.   Histogram  of  first  order  Markov  residuals 
from  logged  rainfall  anomalies  of  data  set  SC 
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This  data  set  yielded  a  chi-square  value  of  273.95 
for  17  degrees  of  freedom.   This  is  equivalent  to  a  signif- 
icance of  zero  plus. 
2.   Winter  Series 
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Figure  75a.   Years  1  -  12  of  winter  months  only, 
logged  rainfall  anomalies  of  data  set  SC. 
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Figure  75b.  Years  13  -  36  of  winter  months  only, 
logged  rainfall  anomalies  of  data  set  SC 
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Figure  75c.   Years  37  -  4  8  of  winter  months 
only,  logged  rainfall  anomalies  of  data  set  SC 
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Figure  76.   Correlogram  of  winter  months  only, 
logged  rainfall  anomalies  from  data  set  SC 
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Figure  77.   Partial  correlogram  of  winter  months 
only,  logged  rainfall  anomalies  from  data  set  SC 
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Figure  78a.   Years  -  24  first  order  Markov  residuals 
of  logged  rainfall  anomalies ,  for  winter  months  only, 
data  set  SC 
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Figure  78b.   Years  25  -  48,  first  order  Markov 
residuals  of  logged  rainfall  anomalies  for  winter 
months  only,  data  set  SC 
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Figure  79.   Correlogram  of  first  order  Markov  residuals 
of  logged  rainfall  anomalies  from  winter  months  only, 
data  set  SC. 
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Figure  80.   Lag  one  plot  of  first  order  Markov  residuals 
from  logged  rainfall  anomalies  of  winter  months  only, 
data  set  SC 
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Figure  81.   First  order  Markov  residuals  versus  lag 
one  data  point  from  logged  rainfall  anomalies  of 
winter  months  only,  data  set  SC 


TABLE  29 

ACTUAL   AND   EXPECTED   NUMBER   OF    TURNING   POINTS 
AND   ACTUAL   AND   EXPECTED   PHASE    FREQUENCIES    FROM 

THE    FIRST    ORDER  MARKOV    RESIDUALS    OF    THE 
LOGGED    RAINFALL   ANOMALIES    OF    THE    WINTER  MONTHS 
ONLY,    DATA   SET    SC 


NUMBER  OF    TURNING   POINTS    =    32 

E[P]    =  30.667  V[P]    =    15.39 

PHASE    LENGTHS 

D  OBS.  E[*] 

1  21  18.8 

2  8  8.1 

3  3  2.3 

4  0.5 

5  0  0.0 

6  0  0.0 

7  0  0.0 

8  0  0.0 

9  0  0.0 
10  0  0.0 

TOTALS  32  29.7 
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TABLE  30 

GENERAL  STATISTICS  OF  FIRST  ORDER 
MARKOV  RESIDUALS  FROM  LOGGED 
RAINFALL  ANOMALIES  OF  WINTER  MONTHS 
ONLY,  DATA  SET  SC 


Moments 

Mean 

.000 

Variance 

.305 

Skewness 

.015 

Kurtosis 

-.363 

Percentiles 


Minimum 

Lower   Sixteenth 

Lower  Eight 

Lower   Quartile 

Median 

Upper   Quartile 

Upper  Eight 

Upper   Sixteenth 

Maximum 

STANDARDIZED    N< 

-1.513 

-  .872 

-  .663 

-  .359 

-  .026 
.355 
.682 
.882 

1.338 
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Figure  82.   Standardized  normal  plot  of  first  order  Markov 
residuals  from  logged  rainfall  anomalies  of  winter  months 
only,  data  set  SC 
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Figure  83.  Histogram  of  first  order  Markov  residuals 
from  logged  rainfall  anomalies  of  winter  months  only, 
data  set  SC 


The  chi-square  was  calculated  at  16.60  for  17 
degrees  of  freedom.   This  is  significant  at  the  .482 
level  thus  indicating  probable  normality. 
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IV.   VALIDATION  OF  FIRST  ORDER  MARKOV  MODELS 

A.   THEORY 

The  general  model  proposed  by  a  first  order  Markov 
process  is,  as  stated  before; 

Z.  =  PZ   +  afc  IV. 1 

where   {a  }  are  distributed  iid  N(0,  a2).   To  validate 
t  a 

this  model,  preferably  independent  data  should  be  subjected 
to  the  model,  and  an  analysis  of  the  residual,  or  forecast 
errors,  made. 

As  stated  previously,  years  19  75  through  19  80  were 
reserved  for  the  purpose  of  validation.   The  method  of 
validation  was  to  use  the  model  to  construct  a  series  of 
one  step  ahead  forecasts.   Let  e. (1)   be  the  error  in  a 
forecast  of  time   t+1   from  the  model  at  time  t.   Then  the 
minimum  mean  squared  error  forecast  (see  Box  and  Jenkins)  is; 

et(l)  =  Zfc  -  PZt_x  IV. 2 

If  the  model  is  correct,  the  sequence   {e  (1) }   will  be 

independent  normally  distributed  with  mean  zero  and  variance 

cf2.   In  the  following  sections  the  models  are  applied  to  the 
a 

reserved  data  sets  (which  may  also  be  found  in  the  appendixes) , 
and  these  forecast  errors  are  calculated.   The  forecast  errors 
are  then  analyzed  to  determine  if 

(1)  The  errors  are  serially  independent 
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(2)  The  errors  are  distributed  as  normal  random 

variables  with  mean  zero,  and  variance   a2. 

a 

Since  the  residual  analysis  of  the  twelve  month  model 
already  indicates  a  poor  fit,  the  twelve  month  model  will 
not  be  validated.   Only  the  winter  month  models  will  be 
checked  for  validity. 

B.   DATA  SET  RN 

Figures  84  (Raw  data)  and  85  (Logged  anomalies)  display 
the  reserved  data  set.   The  logged  anomalies  were  formed  by 
removing  the  means  of  the  analyzed  data,  Table  4,  not  the 
means  of  the  logged  reserved  data.   This  was  done  to  remove 
any  bias  from  the  validation. 
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Figure  84.   Reserved  rainfall  data  for  data  set  RN 
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Figure  85.   Logged  rainfall  anomalies  of  reserved  data  set  RN 


The  forecast  errors,  Figure  86,  their  correlogram, 
Figure  87,  and  independence  tests,  Table  31  indicate  that 
the  errors  are  indeed  independent. 
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Figure  86 .   Forecast  errors  from  first  order  Markov 
model  applied  to  winter  months  of  reserved  data  set  RN 
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Figure  87.   Correlogram  of  forecast  errors  from  first 
order  Markov  model  applied  to  winter  months  of  reserved 
data  set  RN 
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TABLE  31 

ACTUAL  AND  EXPECTED  NUMBER  OF  TURNING 
POINTS  AND  ACTUAL  AND  EXPECTED  PHASE 

FREQUENCEIS  FOR  THE  FORECAST  ERRORS  OF 
THE  FIRST  ORDER  MARKOV  MODEL  APPLIED 

TO  THE  WINTER  MONTHS  OF  RESERVED  DATA 

SET  RN 

NUMBER  OF  TURNING  POINTS  =  32 
E[P]  =  30.667        V[P]  =  15.39 

PHASE  LENGTHS 


D 

OBS. 

E[*] 

1 

21 

18.8 

2 

8 

8.1 

3 

3 

2.3 

4 

0 

.5 

5 

0 

0.0 

6 

0 

0.0 

7 

0 

0.0 

8 

0 

0.0 

9 

0 

0.0 

10 

0 

Q.O 

TOTALS 

32 

29.7 

The  normality  of  the  forecast  errors  is  addressed  by 
Table  32  (Statistics) ,  Figure  88  (Normal  plot) ,  89  (Histogram) , 
and  a  simple  chi-squared  test.   The  chi-squared  was  calculated 
as  7.82  with  5  degrees  of  freedom  which  is  significant  at  the 
.16  7  level.   However,  the  normality  of  the  errors  is  somewhat 
questionable  due  to  the  other  displays. 
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TABLE  32 

GENERAL  STATISTICS  OF  FORECAST  ERRORS  FROM 
THE  FIRST  ORDER  MARKOV  MODEL  APPLIED  TO  THE 
WINTER  MONTHS  OF  RESERVED  DATA  SET  RN 

Moments 


Mean 
Variance 
Skewness 
Kurtosis 


.001 

.259 

-.276 

-.955 


Percentiles 


Minimum 

Lower  Sixteenth 

Lower  Eight 

Lower  Quartile 

Median 

Upper  Quartile 

Upper  Eight 

Upper  Sixteenth 

Maximum 


■1.156 

•  .791 
■    .634 

•  .397 
.069 
.409 
.569 
.614 
.855 
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Figure  88.  Standardized  normal  plot  of  forecast  errors 
from  the  first  order  Markov  model  applied  to  the  winter 
months  of  reserved  data  set  RN 
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Figure  89.   Histogram  of  forecast  errors  from  the 
first  order  Markov  model  applied  to  the  winter 
months  of  reserved  data  set  RN 
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C.   DATA  SET  FL 

As  before,  the  similarity  of  results  for  the  different 
data  set  allows  the  analysis  to  be  portrayed  using  the 
displays  only. 
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Figure  90.   Reserved  rainfall  data  for  data  set  FL 
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Figure  91.   Logged  rainfall  anomalies  of 
reserved  data  set  FL. 
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Figure  92.   Forecast  errors  from  the  first  order 
Markov  model  applied  to  the  winter  months  of 
reserved  data  set  FL 
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Figure  9  3.   Correlogram  of  forecast  errors  from  first 
order  Markov  model  applied  to  the  winter  months  of 
reserved  data  set  FL 

TABLE  33 

ACTUAL   AND    EXPECTED   NUMBER   OF    TURNING    POINTS    AND   ACTUAL 
AND   EXPECTED    PHASE    FREQUENCIES    FROM    THE    FORECAST    ERRORS 
OF    THE    FIRST    ORDER  MARKOV   MODEL   APPLIED    TO    THE    WINTER 
MONTHS    OF    RESERVED    DATA   SET    FL 

NUMBER  OF    TURNING   POINTS    =    34 


E[P]    =    30.667 


V[P]    =    15.396 


PHASE 

LENGTHS 

D 

OBS. 

1 

24 

2 

8 

3 

2 

4 

0 

5 

0 

6 

0 

7 

0 

8 

0 

9 

0 

10 

0 

TOTALS 

34 

E[*] 

18.8 

8.1 

2.3 

.5 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

29.7 
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TABLE  34 

GENERAL  STATISTICS  OF  FORECAST  ERRORS 
FROM  THE  FIRST  ORDER  MARKOV  MODEL  APPLIED 
TO  THE  WINTER  MONTHS  OF  RESERVED  DATA  SET  FL 

Moments 


Mean 

Variance 
Skewness 
Kurtosis 


-.015 

.215 

-.155 

-.970 


Percentiles 


Minimum 

Lower  Sixteenth 

Lower  Eight 

Lower  Quartile 

Median 

Upper  Quartile 

Upper  Eight 

Upper  Sixteenth 

Maximum 


•1.031 

■  .724 
•    .549 

■  .396 
.049 
.369 
.479 
.531 
.904 
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Figure  94.  Standardized  normal  plot  of  forecast 
errors  from  the  first  order  Markov  model  applied 
to  the  winter  months  of  reserved  data  set  FL 
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Figure  95.   Histogram  of  froecast  errors  from  the 
first  order  Markov  model  applied  to  the  winter 
months  of  reserved  data  set  FL 


The  chi-squared  statistic  was  calculated  as  12.5  8 
with  7  degrees  of  freedom,  thus  yielding  a  significance 
level  of  0.083.   This  statistic  and  the  displays  imply 
that  the  data  are  only  marginally  normal  if  at  all. 
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D.   DATA  SET  SC 
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Figure  96.   Reserved  rainfall  for  data  set  SC 
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Figure  97.   Logged  anomalies  of  reserved  data  set  SC 
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Figure  98.  Forecast  errors  from  first  order 
Markov  model  applied  to  the  winter  months  of 
reserved  data  set  SC 
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Figure  99.   Correlogram  of  forecast  errors 
from  first  order  Markov  model  applied  to  the 
winter  months  of  reserved  data  set  SC 
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TABLE    35 

ACTUAL    AND   EXPECTED   NUMBER   OF    TURNING   POINTS    AND 
ACTUAL   AND   EXPECTED    PHASE    FREQUENCES    FROM   THE 
FORECAST    ERRORS    OF    THE    FIRST    ORDER  MARKOV   MODEL 
APPLIED    TO    THE    WINTER   MONTHS    OF    RESERVED    DATA    SET    SC 


NUMBER   OF    TURNING   POINTS 
E[P]    =    30.667 
PHASE    LENGTHS 


=    34 


V[P]    =    15.396 


D 
1 
2 
3 
4 
5 
6 
7 
TOTALS 


OBS. 

El*] 

24 

18.8 

8 

8.1 

2 

2.3 

0 

.5 

0 

0.0 

0 

0.0 

0 

0.0 

34 

29.7 

TABLE    36 

GENERAL    STATISTICS    OF    FORECAST    ERRORS 
FROM   THE    FIRST    ORDER   MARKOV   MODEL   APPLIED 
TO    THE    WINTER  MONTHS    OF    RESERVED    DATA   SET    SC 

Moments 


Mean 

Variance 
Skewness 
Kurtosis 


-.003 

.296 

-.298 

-.413 


Percentiles 


Minimum 

Lower  Sixteenth 

Lower  Eight 

Lower  Quartile 

Median 

Upper  Quartile 

Upper  Eight 

Upper  Sixteenth 

Maximum 


■1.386 

■  .846 

■  .559 
•    .434 

.031 
.414 
.588 
.735 
1.005 
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Figure  100.   Standardized  normal  plot  of  forecast 
errors  from  the  first  order  Markov  model  applied  to 
the  winter  months  of  reserved  data  set  SC 
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Figure  101.   Histogram  of  forecast  errors  from  the 
first  order  Markov  model  applied  to  the  winter 
months  of  reserved  data  set  SC 
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The  chi-squared  statistic  was  calculated  as  9.71  with. 
7  degrees  of  freedom,  thus  yielding  a  significance  level 
of  .205. 

E.   CONCLUSIONS 

The  application  of  a  Markovian  model  was  indicated  by 
the  apparent  dependence  of  adjacent  months  and  the  apparent 
lack  of  dependence  at  any  other  lag.   The  preceeding  sub- 
sections, however,  indicate  that  the  first  order  Markovian 
model  is  weak  at  best. 

The  structure  of  the  data,  visually,  still  points 
toward  some  sort  of  underlying  order.   The  following 
sections  attempt  to  discover  this  order. 
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V.   2x2  TABLES 

A.   THEORY 

As  seen  in  sections  III  and  IV,  the  classical  ARMA  time 
series  approach  does  not  seem  to  adequately  describe  the 
data.   Another  technique  used  to  explore  possible  relation- 
ships is  the  2x2  contingency  tables. 

The  idea  to  be  explored  is  whether  or  not  some  subset  of 
the  data,  to  be  called  the  control,  may  be  used  to  predict 
in  some  way  the  behavior  of  another  subset  of  the  data,  to 
be  called  the  complement.   Here,  the  data  are  reduced  from 
monthly  observations  to  yearly  observations  as  described 
below. 

Let   X  be  the  subset  of  a  year,  to  be  called  the  control, 
and  let  Y  be  the  subset  to  be  called  the  complement.   It  is 
necessary  that  XQY  =  0 ;  that  is,  the  intersection  of  these 
two  sets  is  empty.   The  data  are  then  compared  for  some 
quality  in   X   and  for  some  quality  in   Y.   The  question  is 
then:   does  the  presence  (or  absence)  of  the  quality  in   X 
affect  the  presence  (or  absence)  of  the  quality  in   Y?   An 
example  of  a  typical  table,  is  shown  below  in  Figure  102. 
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Figure  102.   Typical  2x2  contingency  table 


The  table  elements,   n.  .,   represents  the  number  of  years 
which  display  quality   i   in  the  control  and  quality   j   in 
the  complement.   The  marginal  entries   n..   and  n..   represents 
the  numbers  of  years  for  which  the  control  has  quality   i  and 
the  number  of  years  the  complement  has  quality   j   respectively, 
The  overall  number  of  years,   n..,    is  in  the  lower  right  of 
the  table. 

Brownlee  [Ref.  6]  contains  a  very  good  discussion  of  the 
theory  and  use  of  2x2  contingency  tables.   Using  the  notation 
of  Brownlee,  let   9. .   be  the  probability  that  any  given  year 
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will  have  a  control  quality   i   and  a  complement  quality   j . 
Then  estimates  of  the   0 . . *s   are 

/N. 

0 .  .  =  n .  ./n. . 
id     2.y 

8.   =  n.  /n. .  V.l 

9.j  =  n../n 

If  the  control  and  complement  are  independent, 

e..=e..e..  V.2 

id    i   : 

These  simple  assumptions  allow  for  a  thorough  investigation 
of  the  possible  interrelationships  within  the  data  sets. 

Another  way  to  view  the  assumption  of  independence  is 
through  the  use  of  proportions.   Thus,  if  the  basic  division 
is  made  via  the  quality  of  the  control,  the  the  proportions 
P,  =  n^/n^.    (respectively   P2  =  n21/n2.)       V.3 

represent,  in  words,  the  proportion  of  the  years  that  have 
quality   1   in  the  control  and  have  quality   1   (respectively 
2)  in  the  complement. 

The  question  of  independence  may  be  approached  in  several 
ways  as  described  below. 

1.   Fishers  Exact  Test 

A  test  for  the  significance  of  any  dependence  was 
proposed  by  Fisher  in  the  case  in  which  the  marginal  totals, 
n..,n..,  and  n. .   are  known  a  priori  (cf.  Brownlee  [Ref.  6]). 
To  draw  from  Brownlee,  knowledge  of  the  marginals  and  n., 
gives  knowledge  of  all  the  other  elements  of  the  table.   The 
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probability  of  the  event  of  having  exactly  n.,   years  that 
display  quality  1  in  both  the  complement  and  the  control  is; 


nm      -«  \    -      ni     In.  ,  In.  0  in0.  ! 

P(N11-n11)  -  _JL. 1 2_J v>4 

11.  .  •"11-"12  "21   22" 


A  test  may  then  be  applied,  using  V.4,  to  determine  the 
significance  of  any  dependence.   This  test  is  usually  applied 
by  simply  summing  these  probabilities  in  the  tail  of  the 
distribution  (V.4)  in  the  same  direction  as  the  noted  extreme. 

The  usual  procedure  to  provide  a  two-sided  test  of 
significance  is  to  double  a  one  sided  figure.   This  procedure 
is  acceptable  due  to  the  symmetric  appearance  of  the  distri- 
bution. 

Under  the  assumptions  of  independence 

n  ,  .n. 


[„        n,.n.-n0.n. 
Nn   ■  _J: Li 
XXJ      n..  (n..-l) 


V.5 


V  Nni    =  ~j_JL2—l  V*6 

L  11J      n..  (n..-l) 

and   the   random  variable      U      defined  as 

■  u    -E[N      ] 


U  = 


>/v[Nn] 


V.7 


is  asymptotically  distributed  as  a  normal  random  variable 
with  mean  zero  and  variance  one.  This  asymptotic  result 
combined  with  a  continuity  correction  yields  a  test  statistic 
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of 


{  |.n,  .  n-9  ~ni  9  n9i  |.-n.  .  /2  }  / n .  . 
U'   =   ~ — — — V.8 


V  n-,  .n.  ,n_  .n.  _ 

The  statistic   U1   may  then  be  used  as  a  test,  using 
standard  normal  tables,  of  the  significance  of  any  variation 

from  the  assumption  of  independence.   It  should  be  noted  at 

2 
this  point  that  if  the  random  variable   U*   is  squared,  U1 

will  be  distributed  as  a  chi-square  with  one  degree  of  freedom. 

The  squaring  of   U*   with  simplifying  algebra  yields  the  Yates 

correction  to  a  standard  chi-squared  goodness  of  fit  statistic 

2    m    {inll  n22  "n12  n21  l-n--/2>    n..       y>9 
(n11+n12)  Cn11+n21) (n12+n22) (n21+n22) 

see  Dixon  and  Massey  [Ref.  5].   This  allows  the  use  of  the 

chi-square  tables  as  an  equivalent  test  to  that  of  V.8. 

2.   Odds 

Subsection  V.A.I  above  deals  with  the  significance 

of  any  observed  interdepedence  between  the  control  and  the 

complement.   The  question  of  the  degree  of  dependence  should 

also  be  addressed.   The  measure  to  be  used  is  the  odds  ratio. 

Using  the  notation  of  Fleiss  [Ref.  3],  a  measure  of  seeing 

quality  1  in  the  complement  Y   may  be 

m      P(Y=H  X=l)   . 
1     P(Y=2|  X-l)  '  V#  U 

this  is  then  the  odds  that  quality  1  will  occur  in  the  comple- 
ment given  that  quality  1  is  present  in  the  control.   In  a 
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similar  manner; 

a    -  pfY=1l  x=2)  .. 

U2  P(Y=2 |  X=2)  V-11 

is  the  odds  that  quality  1  will  occur  in  the  complement  given 
that  quality  2  was  observed  in  the  control.   The  currently 
most  often  used  measure  is  the  odds  ratio   oo,   or 

Ql  V.12 

"  =—2 

Note  that,  if  the  appearance  of  quality  1  in  the  complement 
is  independent  of  whether  or  not  it  appears  in  the  control, 
then  a)  =  1.   While   oo>l   implies  that  the  odds  of  the 
complement  having  quality  1,  given  that  quality  1  was  observed 
in  the  control,  are  greater  than  the  odds  of  the  complement 
having  quality  2.   This  would  indicate  that  the  control  would 
be  some  sort  of  predictor  for  the  complement,  relative  to  the 
selected  qualities. 

In  the  same  continuity  correcting  spirit,  as  was  used 
with  the  Yates  chi-square,  an  estimate  for   u>  may  be  obtained 
from  the  Table  as 


a  =  0  =   (nll+'5>  <n22+-5) 


(n12+.5) (n21+.5) 


V.13 


with  a  standard  error  of 


s.e.(0)=0\/ — =r— =■  +  — ±— r  +  — ±-5*.  V.14 

Vnn+.5    n12+.5    n21+<5  +  n22+.5 

The  natural  logarithm  of  this  odds  ratio  will  be  discussed 
more  fully  in  section  VI . 
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B.   ANALYSIS 

The  theory  of  subsection  A  above  is  applied  to  the  three 
data  sets  as  discussed  below.   The  control  is  typically  taken 
as  a  monthly  anomaly,  say  October.   Here,  the  quality  is  taken 
as  either  a  positive  or  a  negative  anomaly.   Thus   X=l   occurs 
when  the  month  of  October  falls  below  its  mean  and   X=2   occurs 
when  it  falls  above  its  mean.   The  complement  consists  of  the 
sum  of  the  rainfall  for  the  succeeding  eleven  months,  or  in 
symbols : 

Xt  =  Rt,l  "  V 


12  ,   N    12 


Yt  = 


=  I       K         -  h-     l    (  I   R*.  )  •  V.15 

L~         t,m    N  .L,         L-  t,m 
m=2     '        t=l   m=2   ' 


Where  it  is  understood  that   X=l  when   X  <  0,   X=2   when 

X,  >  0   and  similarly  for   Y  . 

Various  control  subsets  are  used;  October  through 

September  were  investigated  by  themselves  as  were  all 

adjacent  pairs,  triples,  and  four-tuples  of  months.   For 

an  example,  consider  the  spring  (April,  May,  and  June)  and 

its  complement  (July  through  March) .   In  this  case 

9  N-l    9 

t  =   S  Rt,m  "  i-  J>.     (  S  R4.   )  v-16 

m=7    '     N  t=l   m=7   t,m 


X 


12  6  N-l  12  6 

Yt  =  JlORt'ra  +  mil    Rt+l-»  '  F1  tllUlO^1"  +  JlRt+1'm)  • 
Equations  V.15  and  V.16  imply  that  the  data  are  always 
analyzed  as  deviations  from  the  arithmetic  mean.   Eowever, 
the  data  are  also  analyzed  as  deviations  from  the  median 
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and  the  lower  quartile.   In  the  tables  to  follow,  'A'  refers 
to  both  control  and  complement  having  the  arithmetic  mean 
removed,  'M'  refers  to  both,  control  and  complement  having 
their  respective  medians  removed,  and  ' QL'  refers  to  the 
control  having  the  lower  quartile  removed  while  the  median 
was  removed  from  the  complement. 

The  first  four  Tables  (.37  through  40)  ,  give  the  signifi- 
cance levels  of  observed  departures  from  independence  of  the 
control  and  complement.   Only  those  values  having  a  Yates 
corrected  chi-square  of  greater  than  1.00  are  listed.   The 
entries  represent  the  two-tailed  probability  of  a  random 
deviation  in  excess  of  that  observed.   Although  the  cut  off 
criterion  was  the  Yates  chi-square,  the  agreement  between 
its  probability  and  that  obtained  from  the  Fisher  exact  and 
normal  tests  did  not  differ  in  the  first  two  decimal  places. 

TABLE  37 

SIGNIFICANCE  OF  OBSERVED  DEPARTURES  FROM  INDEPENDENCE  OF 
SINGLE  MONTH  CONTROL  VERSUS  SUCCEEDING  ELEVEN 
MONTH  COMPLEMENTS 

Data  set  RN        FL  SC 

Differentiator   A   M   QL   A   M   QL   A   M   QL 

Control 

October  .18   .22  .24  .31       .15 

November  . 10  . 12  .28 

December 

January  .14   .30  .19  .002  .04  .12  .21   .19 

February  .  11 

March.  .13 

April  .18 

May  . 14 

June  .24 

July  .31  .21 

August 

September  .14  .12 
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TABLE  38 

SIGNIFICANCE  OF  OBSERVED  DEPARTURES  FROM  INDEPENDENCE 
OF  PAIRS  OF  MONTHS  VERSUS  SUCCEEDING 
TEN  MONTH  COMPLEMENTS 


Data  Set 

RN 

FL 

SC 

Differentiators 

A    M 

QL 

A 

M 

QL 

A 

M 

QL 

Control 

Oct+Nov 

.18 

.09 

.24 

Nov+Dec 

.10 

.12 

.15 

.31 

.08 

Jan+Feb 

.30  .30 

.01 

.02 

.12 

.11 

.19 

Feb+Mar 

.11 

Mar+Apr 
Apr+May 
May+Jun 
Jun+Jul 

.14 

.12 

.24 

Jul+Aug 
Aug+Sep 

.19 

.18 

TABLE    39 

SIGNIFICANCE    OF    OBSERVED    DEPARTURES    FROM    INDEPENDENCE 
OF    TRIPLES    OF   MONTHS    VERSUS    SUCCEEDING 
NINE-MONTH    COMPLEMENTS 

SC 

M         QL         A         M         QL 


19 


Data  Set 

RN 

FL 

Differentiators 

A    M 

QL    A 

Control 

Oct+Nov+Dec 

Nov+Dec+Jan 

Dec+ Jan+Feb 

~ 

.31 

Jan+Feb+Mar 

Feb+Mar+Apr 

Mar+ Apr+May 

Apr+May+Jun 

.30 

May+Jun+Jul 

.14 

.30 

Jun+Jul+Aug 

Jul+Aug+Sep 

.31    .28 
22 


.18 
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TABLE    40 

SIGNIFICANCE    OF    OBSERVED    DEPARTURES    FROM 

INDEPENDENCE    OF    FOUR-TUPLES    OF   MONTHS    VERSUS 

SUCCEEDING   EIGHT-MONTH    COMPLEMENTS 

Data    set  RN  FL  SC 

Differentiator      A        M        QL        A        M        QL        A        M        QL 

Control 

Oct+Nov+Dec+Jan  .26 

Nov+Dec+Jan+Feb      .22 

Dec+Jan+Feb+Mar      .22 

Jan+Feb+Mar+Apr  .31 

Feb+Mar+Apr+May   .22  .27 

Mar+Apr+May+Jun  .20  .31 

Apr+May+Jun+Jul  .30  .28 

May+Jun+Jul+Aug       . 14 

Jun+Jul+Aug+Sep 


Several  choices  for  predictors  are  suggested  in  the 
previous  tables.   However,  the  apparent  strongest  candidate 
for  a  predictor  is  January.   The  control  of  January  by 
itself  and  January  paired  with  December,  are  the  most 
consistently  significant  entries.   Tables  41  below  gives 
the  odds  ratio,  V.13,  for  January,  and  January  and  December, 
as  controls. 
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TABLE  41 

ODDS  RATIO  OF  JANUARY  VERSUS  FEBRUARY 

THROUGH  DECEMBER  AND  JANUARY  PLUS 

DECEMBER  VERSUS  FEBRUARY  THROUGH  NOVEMBER 


Differentiator  A  M  QL 
Data  set  RN 

January  4.59  3.15  5.13 

Jan+Dec  3.15  3.15  2.01 

Data  set  FL 

January  10.71  4.72  4.30 

Jan+Dec  7.42  6.02  4.30 

Data   set  SC 

January  2.45  2.49  2.23 

Jan+Dec  3.00  2.49  1.44 


At  this  point  in  the  analysis  it  was  decided  to  explore 
more  fully  the  power  of  January  as  a  predictor.   It  should 
be  stated  that  other  possibilities  for  predictors  are 
suggested  by  the  tables,  but  time  did  not  allow  an  exhaustive 
study  of  all  of  these. 

C.   OTHER  RESULTS 

The  results  of  section  V.B  suggest  that  a  more  detailed 
analysis  of  January  as  a  predictor  is  in  order.   The  first 
method  tried  for  this  was  ordinary  least  squares  regression 
of  the  rainfall  total  in  January  versus  the  total  for  February 
through  December.   This  is  the  model  below. 

Let  Xt  -  Rt/4 


12 

Y.  =   J  R.    +   y  RL1 ,  V.17 

,m 


J.4  3 

'  =  Jr.  +  y  r.  . , 

t     **_  t,m     S  t+1,: 
m=5       m=l 
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then  assume  that 


Y.  =   a  +  ex.  +  e. 
t  t    t 


V.18 


as  the  standard,  linear  model  where  tex.}    are  assumed  to  be 

independent  and  identically  distributed  with  mean  zero  and 

variance  ae2.   If  the  predictability  of  January  is  strong, 

this  model,  V.18,  may  result  in  a  good  fit  of  the  data. 

Table  42  below  is  the  resulting  ANOVA  for  this  regression. 

As  may  easily  be  seen,  the  model  does  not  appear  to  have  any 

significance. 

TABLE  4  2 

ANOVA  FOR  REGRESSION  OF  SIMPLE 
LINEAR  MODEL  FOR  ALL  DATA  SETS 


R-squared  =  .1033 

Standard  error  of  estimate  =  4.1574 


RN 


SOURCE 

Total 

Regression 
Jan-Control 
Residual 

Variable 
Alpha 
Beta 


Alpha 
Beta 


DF 

22 

1 

1 

21 

Coefficient 
.2178 
.5993 


SS 

404.798 
41.821 
41.821 

362.977 


MS 


41.821 
41.821 
17.285 


2.42 
2.42 


Standard  error   T 
.878         .25 
.385        1.56 


90%  Confidence  Limits 

Lower  limit       Upper  limit 
-1.1817  1.6174 

-  .0148  1.2133 
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R-squared  =  .2066 

Standard  error  of  estimate  =  4.28Q4 


FL 


AOV 


SOURCE 

DF 

SS 

M£ 

■ 

F 

i 

Total 

35 

785. 

196 

Regression 

1 

162. 

,252 

162. 

252 

8. 

86 

Jan -Control 

1 

16  2. 

252 

162. 

252 

8. 

86 

Residual 

34 

622. 

,945 

18. 

322 

Variable 

Coefficient 

Standard 

error 

T 

Alpha 

-.2670 

a 

715 

-, 

,37 

Beta 

1.0095 

i 

339 

2. 

,98 

Alpha 
Beta 


90%  Confidence  Limits 
Lower  limit       Upper  limit 
-  .2670  .8543 

1.0095  1.5417 


R-squared  =  .0426 

Standard  error  of  estimate  =  6.150  8 


SC 


SOURCE 

Total 

Regression 
Jan-Control 
Residual 

Variable 

Alpha 
Beta 


Alpha 
Beta 


DF 

46 

1 

1 
45 

Coefficient 

-.2771 
.4355 


AOV 
SS 

1778.131 
75.687 
75.687 

1702.444 


MS 


75.687 
75.687 
37.832 


Standard  error 

.898 
.308 


2.00 
2.00 


-.31 
1.41 


9  0%  Confidence  Limits 
Lower  limit      Upper  limit 
-1.6766  1.1224 

-  .0446  .9157 


The  same  model  with  means  removed,  V.19,  below  was  tried 

and,  although  slightly  better,  is  still  not  strong. 

Y  -Y.  =   a  +  3(X.-X. )  +  e^.  V.19 

t  t        t 
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TABLE    4  3 

ANOVA   FOR    REGRESSION    OF    SIMPLE    LINEAR 
MODEL   WITH   MEANS    REMOVED    FOR   ALL    DATA   SETS 


R-squared  =    .1033 

Standard  error   of   estimate   =   4.1575 


AOV 


RN 


SOURCE 


DF 


SS 


MS 


Total  22 

Regression  1 

Jan-Control  1 

Residual  21 


404.798 
41.821 
41.821 

362.977 


41.821 
41.821 
17.285 


2.42 
2.42 


Variable 
Alpha 
Beta 


Coefficient 
11.7215 
.5993 


Standard  error  T 

1.851  6.33 

.385  1.56 


Alpha 


90%  Confidence  Limits 
Lower  limit         Upper  limit 
8.7236  13.2483 

-.0148  1.2133 


R-squared  =  .174  7 

Standard  error  of  estimate  =  4.4276 


FL 


AOV 


SOURCE 


DF 


SS 


MS 


Total  35 

Regression  1 

Jan-Control  1 

Residual  34 


807.618 
141.108 
141.108 
666.510 


141.108 

141.108 

19.603 


7.20 
7.20 


Variable 
Alpha 
Beta 


Coefficient 
10.9860 
.9414 


Standard  error      T 
1.442  7.62 

.351  2.68 


Alpha 
Beta 


90%  Confidence 
Lower  limit 
8.7236 
.3909 


Limits 

Upper  limit 
13.2483 
1.4920 
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R-squared  =  .5617 

Standard  error  of  estimate  =  .7589 


SC 


SOURCE 


Alpha 
Beta 


DF 


AOV 


ss 


MS 


Total 

Regression 
Jan-Control 
Residual 

7 
1 
1 
6 

7.884 
4.428 
4.428 
3.456 

4.428      7.69 
4.428      7.69 
.576 

Variable 
Alpha 
Beta 

Coefficient 
-  .1946 
1.7710 

Standard  error       T 
.287           -.68 
.637           2.77 

90%  Confidence  Limits 
Lower  limit        Upper  limit 
-.7040  .3148 

.6362  2.9058 


The  amount  of  rainfall  in  January  does  not  appear  to  be 
a  strong  predictor  for  the  amount  of  rainfall  in  February 
through  December.   This  seems  to  indicate  that  the  relation- 
ships between  January  rainfall  and  rainfall  during  the  next 
eleven  months  is  not  as  strong  as  expected.   However,  a 
further  technique  is  available,  that  of  log-odds  and  logistic 
regression,  which  are  the  subjects  of  the  next  section. 
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VI.   LOGISTIC  ANALYSIS 

A.   THEORY 

The  logistic  analysis  to  be  described  in  the  section  was 
developed  from  Gaver  [Ref.  2]  and  Fleiss  [Ref.  3].   This 
analysis  derives  from  the  model  of  1.5  as  stated  in  the 
introduction. 

The  basic  approach  is  to  view  the  complement  as  having  a 
binary  representation,  with  success  being  defined  as  a  comple- 
ment above  its  mean  (see  equations  1.3  and  1.4)  and  failure 
as  the  complement  below  its  mean.   The  problem  then  is  to  find 
the  conditional  probability  of  a  success  (the  complement  being 
above  its  mean  for  a  year)  given  that  the  control  (January 
rainfall)  takes  on  a  particular  value. 

The  control  is  now  taken  to  be  the  logged  rainfall  anomaly 
of  January  and  is  found  in  equation  1.2  repeated  below. 

xt  -  &(Rt,4'  -  irJ^V^  VI-X 


If  the  probability  of  success,  given  X   is  written  as 

P (Success | X  )  -  9  VI. 2 

a  superficially  attractive  model  for  9    is 

9t  =   a  +6  Xt  +  efc  VI. 3 

This  model  has  two  difficulties,  the  worst  of  which,  is 
that  probabilities  of  greater  than  one  or  less  that  zero  are 
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allowed.   Secondly,  the   6    are  not  available  in  proportion 

form  to  fit  the  model. 

Initially,  the  problem  of  estimating   9.   is  approached 

by  grouping  the  data.   Section  II  indicated  that  each  year 

seemed  to  be  independent  of  the  next  and  that  there  was  no 

trend.   This  allows  for  the  ordering  of  the  X,_ '  s  into  their 

3         t 

order  statistics   X ,  .  's.   Once  the  ordering  has  been  done, 
non  overlapping  groups  of  arbitrary  size  may  be  formed  as 
shown.   Let  X, ,X~ /Xo* • • • >X, ?   ^e  a  series  °f  12  years  from  an 
arbitrary  data  set,  with  associated  order  astatistics 

X(1)'X(21'X(3)',,,'X(12) '   T^e  "^  groups  of  size  three  are 
desired,  the  data  are  partitioned  below. 


X 


(1) '  X(2) 'X(3) lX(4) /X(5)  'X(6) I *• * lX(10) /X(11) ,X(12) 


Given  these  groups,  let  X.   be  a  measure  of  location  for 

the  j    group.   This  analysis  used  the  median,  therefore 

X.  =  X,-.,,*.   Also  associated  with  each  group  is  R.  (not 

rainfall) ,  the  number  of  success  in  group   j ,   and  n . ,   the 

number  of  elements  in  group  j.   From  this  set  up,  the  required 

probabilities  may  be  estimated  as; 

9 .  =  R . /n .  VI . 4 

3  3      3 

A  solution  to  the  first  problem,  that  of  the  model  yield- 
ing probabilities  outside  of  (0,1) ,  is  to  use  the  log  odds, 
instead  of   9 .   where; 

Log  odds  =  <J>  .=£n  I  ^q—  j  VI. 5 
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which-  is  equivalent  to  the  logarithm  of  the  odds  ratio  as 
given  in  V.13.   Gaver  [Ref.  2 J  suggests  that  a  correction  of 
.5  be  applied  to  guard  against  the  problem  of  0  and  1  within 
the  logarithm  and  to  reduce  the  bias.   The  statistic  then 

becomes; 

/ ft  +. «5  \ 

VI.  6 


The  temptation  is  to  go  directly  to  the  model 

$j  =  a  +  3    X. ;  VI. 7 

yet  V[9.   =  9.(l-9.)/n.   which  is  not  constant.   This 
suggests  the  need  for  a  weighting  scheme. 

The  weighting  scheme  used  was  that  of  iteratively 
reweighted  least  squares,  (IRWLS) ,  using  the  bi-weights. 
This  method  is  dicussed  in  detail  in  Mosteller  and  Tukey 
[Ref.  13]. 

Although  grouping  of  the  data  and  the  model  of  VI. 7 
provide  adequate  representation  of  the  underlying  structure, 
the  logistic  model  itself,  1.5,  when  viewed  through  the  eyes 
of  maximum  likelihood  theory  can  yield  more  insight. 

The  model  is  assumed  to  be 

ea+6  Xt 


1+  ea"BXt  VI'8 


where  the   X   are  independent. 
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The  likelihood  function  is  then 


N 

YL  /  \     1-Y 


L(X,Y;a,B)    =    f]  ea+BXt  t 


t=l        \l  +   ea+eXt/  \l   +   ea+BXt 


eaXt+6XtYt 


n 


N 
t=l      n    (l+ea+6Xt) 

t=l 


and   the   log-likelihood   is 

N  N 

L(X,Y;a,3)    =    a    J   Y      +         $£      X.Y, 
t=l    C  t=l      r    C 


t=l    t       t       t=l    t    \l    +    tyj 


and 


t=iTi~+ljj^)2  t=i  (i  +  i>t)2 


xt  »t       -!  xt  »t . 

t=l    (1+    ^t)2  t=l    (l+*t) 

where  ^0    „ 

ip      =   e  t. 


t 


V.9 


8       ,       ,,  a+BX.,  V.10 

-2,      £n    (1  +   e  t) 

t=l 
The    gradient,    and  Hessian   of  VI . 9    are 


aa     ^    t     ti1   Vr+-T£. 


V.ll 
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A  simple  way  to  solve  for   a  and   3   is  to  use  Newtons 
Method  as; 

(a)  =  (a)         -    H71  V  fcL 

\3/k+l    W  k 

since  all  necessary  elements  may  be  calculated  in  one  pass  of 
the  computor  algorithm. 

One  beneficial  byproduct  of  the  maximum  likelihood 
approach  is  the  asymptotic  information  matrix,   H,  .   Gaver 
[Ref.  2]  states  that  the  diagonal  elements  of  this  matrix 
provide  good  estimates  of   V[a]  and  V[3]   under  assumptions 
of  normality. 

B.   ANALYSIS 

1.   Grouped  Data 

The  first  approach  taken  was  that  of  grouping  the 
data  as  described  above.   Groups  of  3,  4,  and  5  were  used, 
as  were  two  separate  methods  of  regression,  ordinary  least 
squares  (OLS )  ,  and  iteratively  reweighted  least  squares 
(IRWLS) .   Tables  44  (RN) ,  45  (FL) ,  and  46  (SC)  present  the 
data  and  Tables  47  (OLS) ,  and  48  (IRWLS)  present  the  results 
of  the  regressions. 
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TABLE  4  4 

DATA  SET  RN,  LOGGED  JANUARY  ANOMALIES  AND 
SUCCESSES  FOR  GROUPED  AND  UNGROUPED  FORMS 


UNGROUPED 

t 

X        Y 

1 

-1.22      0 

2 

-1.09      0 

3 

-  .66      0 

4 

-  .58      0 

5 

-  .46      0 

6 

-  .36      1 

7 

-  .31      0 

8 

-  .18      1 

9 

-  .17      0 

10 

-  .05      1 

11 

.01      1 

12 

.01      0 

13 

.15      1 

14 

.16      0 

15 

.24      1 

16 

.28      0 

17 

.36      1 

18 

.46      1 

19 

.48      1 

20 

.50      1 

21 

.51      0 

22 

.94      1 

23 

1.01 

GROUPED 

J 

X 
group=3 

Y 

1 

-1.09 

0 

2 

-  .46 

1 

3 

-  .19 

1 

4 

.01 

2 

5 

.16 

2 

6 

.36 

2 

7 

.50 

2 

8 

.98 
group=4 

2 

1 

-.38 

0 

2 

-.34 

2 

3 

-.03 

2 

4 

.20 

2 

5 

.47 

4 

6 

.94 
group=5 

2 

1 

-.66 

0 

2 

-.19 

3 

3 

.15 

3 

4 

.46 

4 

5 

.94 

2 
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TABLE  45 

DATA  SET  FL,  LOGGED  JANUARY  ANOMALIES  AND 
SUCCESSES  FOR  GROUPED  AND  UNGROUPED  DATA 


UNGROUPED 

GROUPED 

t 

X 

Y 

j 

X 

Y 

1 

-3.41 

1 

group=3 

2 

-1.82 

0 

1 

-1.82 

1 

3 

-1.14 

0 

2 

-1.00 

0 

4 

-1.02 

0 

3 

-  .46 

0 

5 

-1.00 

0 

4 

-  .20 

2 

6 

-  .92 

0 

5 

.11 

1 

7 

-  .47 

0 

6 

.14 

1 

8 

-  .46 

0 

7 

.20 

0 

9 

-  .39 

0 

8 

.34 

2 

10 

-  .28 

1 

9 

.50 

2 

11 

-  .20 

0 

10 

.61 

2 

12 

-  .14 

1 

11 

.75 

3 

13 

-  .04 

1 

12 

1.14 

3 

14 

.11 

0 

group=4 

15 

.11 

0 

1 

-1.48 

1 

16 

.13 

0 

2 

-  .70 

0 

17 

.14 

0 

3 

-  .24 

2 

18 

.16 

1 

4 

.11 

1 

19 

.16 

0 

5 

.16 

1 

20 

.20 

0 

6 

.32 

2 

21 

.25 

0 

7 

.50 

2 

22 

.30 

1 

8 

.71 

4 

23 

.34 

1 

9 

1.11 

4 

24 

.38 

0 

group=5 

25 

.46 

1 

1 

-1.14 

1 

26 

.50 

0 

2 

-  .46 

1 

27 

.51 

1 

3 

-  .04 

2 

28 

.57 

0 

4 

.16 

1 

29 

.61 

1 

5 

.34 

3 

30 

.71 

1 

6 

.57 

3 

31 

.71 

1 

7 

.95 

6 

32 

.75 

1 

33 

.82 

1 

34 

1.08 

1 

35 

1.14 

1 

36 

1.31 

1 
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TABLE  46 

DATA  SET  SC,  LOGGED  JANUARY  ANOMALIES  AND 
SUCCESSES  FOR  GROUPED  AND  UNGROUPED  FORMS 


UNGROUPED 

GROUPED 

t 

X 

Y 

j 

X 

Y 

1 

-4.2  4 

1 

group=3 

2 

-1.41 

0 

1 

-1.41 

1 

3 

-1.38 

0 

2 

-1.07 

0 

4 

-1.08 

0 

3 

-    .96 

1 

5 

-1.07 

0 

4 

-    .55 

1 

6 

-1.06 

0 

5 

-    .26 

2 

7 

-1.02 

0 

6 

-    .10 

1 

8 

-    .96 

1 

7 

.02 

2 

9 

-    .65 

0 

8 

.20 

0 

10 

-    .60 

0 

9 

.28 

2 

11 

-    .55 

1 

10 

.40 

2 

12 

-    .36 

0 

11 

.46 

1 

13 

-    .27 

1 

12 

.56 

2 

14 

-    .26 

0 

13 

.69 

0 

15 

-    .21 

1 

14 

.78 

2 

16 

-    .17 

1 

15 

.97 

1 

17 

-    .10 

0 

16 

1.25 

2 

18 

-    .03 

0 

group=4 

19 

-    .01 

1 

1 

-1.39 

1 

20 

.02 

0 

2 

-1.04 

1 

21 

.12 

1 

3 

-    .58 

1 

22 

.14 

0 

4 

.23 

3 

23 

.20 

0 

5 

-    .02 

1 

24 

.22 

0 

6 

.17 

1 

25 

.23 

0 

7 

.28 

3 

26 

.28 

1 

8 

.42 

2 

27 

.28 

1 

9 

.56 

2 

28 

.40 

1 

10 

.69 

1 

29 

.40 

1 

11 

.78 

1 

30 

.41 

0 

12 

1.03 

3 

31 

.43 

0 

group=5 

32 

.46 

1 

1 

-1.38 

1 

33 

.49 

0 

2 

-    .96 

1 

34 

.56 

1 

3 

-    .27 

3 

35 

.56 

1 

4 

-    .04 

2 

36 

.60 

0 

5 

.20 

1 

37 

.64 

0 

6 

.40 

4 

38 

.69 

0 

7 

.49 

3 

39 

.70 

0 

8 

.69 

1 

40 

.70 

1 

9 

.88 

2 

41 

.78 

1 

10 

1.25 

2 

42 

.81 

0 

43 

.88 

0 

44 

.97 

0 

45 

.97 

1 

46 

1.03 

1 

47 

1.47 

1 
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TABLE  4  7a 

ORDINARY  LEAST  SQUARES  REGRESSION  WITH 
THE  MODEL  OF  EQUATION  VI . 7  FOR  DATA  SET  RN 


R-squared  =  .0423 

Standard  error  of  estimate  =  6.1508 


RN 
GROUP=3 


SOURCE 


DF 


SS 


AOV 


MS 


Total 

Regression 
Jan-Control 
Residual 

46 
1 
1 

45 

1778.131 
75.687 
75.687 

1702.444 

75.687 
75.687 
37.832 

2.00 
2.00 

Variable 
Alpha 
Beta 

Coefficient        Standard  error 
14.5833            1.677 
.4355              .308 

T 
8.69 
1.41 

Alpha 
Beta 


90%  Confidence  Limits 
Lower  limit        Upper  limit 
11.9680  17.1988 

-  .0446  .9157 


R-squared  =  .4068 
Standard  error  of  estimate 


1.2099 


RN 
GROUP=4 


AOV 


SOURCE 


DF 


SS 


MS 


Total 

Regression 
Jan-Control 
Residual 

5 
1 

1 
4 

9.871 
4.016 
4.016 
5.855 

4.016 
4.016 
1.464 

2.74 
2.74 

Variable 
Alpha 
Beta 

Coefficient 
-.1696 
1.7687 

Standard  error 
.517 
1.068 

T 
-.33 
1.66 

Alpha 
Beta 


90%  Confidence  Limits 
Lower  limit        Upper  limit 
-1.1675  .8284 

1.7688  3.8288 
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R-squared  =  .5973 

Standard  error  of  estimate  =  .9995 


RN 
GR0UP=5 


SOURCE 


DF 


AOV 


SS 


MS 


Total  4 

Regression  1 

Jan -Control  1 

Residual  3 


7.442 
4.446 
4.446 
2.997 


4.446 

4.446 

.999 


4.45 
4.45 


Variable 
Alpha 
Beta 


Coefficient 
-.2647 
1.7193 


Standard 
.461 
.815 


error 


T 
-.57 
2.11 


Alpha 
Beta 


90%  Confidence 
Lower  limit 
-1.2358 
.0040 


Limits 

Upper  limit 
.7064 
3.4345 


TABLE  4  7b 

ORDINARY  LEAST  SQUARES  REGRESSION  WITH 
THE  MODEL  OF  EQUATION  VI . 7  FOR  DATA  SET  FL 


R-squared  =  .2862 

Standard  error  of  estimate  = 


.8949 


FL 
GROUP=3 


SOURCE 


DF 


AOV 


SS 


MS 


Total  11 

Regression  1 

Jan-Control  1 

Residual  10 


20.591 
7.998 
7.998 

12.593 


7.998 
7.998 
1.259 


35 
35 


Variable 
Alpha 
Beta 


Coefficient 
-.1459 
1.0521 


Standard  error     T 

.324  -.33 

.417  2.52 


Alpha 
Beta 


90%  Confidence  Limits 
Lower  limit        Upper  limit 
-.6871  .4944 

.3550  1.7492 
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R-squared  =  .5370 

Standard  error  of  estimate  =  1.0484 


FL 
GROUP=4 


AOV 


SOURCE 


DF 


SS 


MS 


Total  8 

Regression  1 

Jan-Control  1 

Residual  7 


16.615 
8.922 
8.922 
7.963 


8.922 
8.922 
1.099 


8.12 
8.12 


Variable 


Coefficient 


Standard  Error 


Alpha 
Beta 


-.1148 
1.3572 


.350 
.476 


-.33 
2.83 


Alpha 
Beta 


90%  Confidence  Limits 
Lower  limit         Upper  limit 
-.7237  .4941 

.5297  2.1847 


R-squared  =  .6  311 

Standard  error  of  estimate  = 


.8809 


FL 
GROUP=5 


AOV 


SOURCE 


DF 


SS 


MS 


Total 

Regression 
Jan-Control 
Residual 

6 

1 
1 
5 

10. 
6, 
6, 
3. 

.519 
.639 
.639 
.880 

6.639 
6.639 

.776 

8.56 
8.56 

Variable 
Alpha 
Beta 

Coefficient 
-.1393 
1.5225 

Standard  Error 
.334 
.521 

T 
-.42 
2.92 

Alpha 
Beta 


90%  Confidence  Limits 
Lower  limit  Upper  limit 

-.7525  .4740 

.5673  2.4777 
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TABLE  4  7c 

ORDINARY  LEAST  SQUARES  REGRESSION  WITH 
THE  MODEL  OF  EQUATION  VI . 7  FOR  DATA  SET  SC 


R-squared  =  .1404 

Standard  error  of  estimate  -  .9901 


SC 
GROUP=3 


AOV 


SOURCE 


DF 


SS 


MS 


Total 

15 

15.965 

Regression 

1 

2.241 

2.241 

2.29 

Jan -Control 

1 

2.241 

2.241 

2.29 

Residual 

14 

13.724 

.980 

Variable 

Coefficient 

Standard 

Error   T 

Alpha 

-.3039 

.249 

-1.22 

Beta 

.5059 

.335 

1.51 

Alpha 
Beta 


90%  Confidence  Limits 
Lower  limit       Upper  limit 
-.7090  .1012 

-.0386  1.0504 


R-squared  =  .1872 

Standard  error  of  estimate  = 


.8959 


SC 
GROUP=4 


AOV 


SOURCE 


DF 


SS 


MS 


Total  11 

Regression  1 

Jan-Control  1 

Residual  10 


9.875 
1.848 
1.848 
8.026 


1.848 

1.848 

.803 


2.30 
2.30 


Variable 
Alpha 
Beta 


Coefficient 
-.2284 
.5447 


Standard  Error       T 
.260  -.88 

.359  1.52 


Alpha 
Beta 


90%  Confidence 
Lower  limit 
-.6618 
.5447 


Limits 

Upper  limit 
.2051 
1.1439 
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R=squared  =  .2862 
Standard  error  of 


estimate  =  . 8949 


SC 
GR0UP=5 


SOURCE 


DF 


ADV 


SS 


MS 


Total  9 

Regression  1 

Jan-Control  1 

Residual  8 


8.976 
2.569 
2.569 
6.407 


569 
569 

801 


3.21 
3.21 


Variable 
Alpha 
Beta 


Coefficient 
-.1838 
.6573 


Standard 
.287 
.367 


Error 


T 
-.64 
1.79 


Alpha 
Beta 


90%  Confidence  Limits 
Lower  limit        Upper  limit 
-.6736  .3061 

.0303  1.2843 


TABLE  4  8 

ITERATIVELY  REWEIGHTED  LEAST  SQUARES 
REGRESSION  USING  BI-WEIGHTS  FOR  THE  MODEL  OF  EQUATION  VI . 7 


GROUP  SIZE  3 


C=9 

a 

Data  sets 

RN 

.048 

FL 

-.063 

SC 

-.116 

C=4 

Data  sets 

RN 

.052 

FL 

-.197 

SC 

.154 

GROUP 

SIZE  4 

C=9 

Data  sets 

RN 

-.103 

FL 

-.093 

SC 

-.169 

GROUP 

SIZE  5 

C=9 

Data  sets 

RN 

-.065 

FL 

-.107 

SC 

-.126 

3 
1.013 

.664 

.351 


1.031 

1.049 

.480 


665 
705 
280 


865 
723 
413 
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2 .   Maximum  Likelihood 

Table  49  displays  the  final  points  for  each  data  set 

along  with  the  inverse  Hessian  at  that  point.   Figures  103 

(RN) ,  104  (FL) ,  and  105  (SC)  are  interesting  in  that  they 

portray  the  contours  of  the  likelihood  functions  for  each 

data  set.   These  contours  show  the  variance  of  the  estimated 

parameters  in  a  graphic  way.   Note  how  data  sets  RN  and  FL 

seem  to  have  some  sort  of  horizontal  ridge  indicating  a 

good  pick  of  the  slope  parameter,  yet  the  contours  of  data 

set  SC  are  almost  circular  about  the  origin  of  the  axes 

indicating  no  significant  difference  from  zero  for  either 

parameter. 

TABLE  49 

MAXIMUM  LIKELIHOOD  ESTIMATES  OF  a  AND  3  ALONG  WITH 
ESTIMATES  OF  THEIR  VARIANCE  FOR  ALL  THREE  DATA  SETS 


Data   set   RN 

/ 

> 

a   =      .062 
3   =   2.918 

"I1      " 

/    .257 
V-.034 

-.043 
1.720 

V[a]  =       .257 

V[0]  =    1.720 

Data    set   FL 
a  =      -.171  , 

3   =         .933 

V[a]  =      .129 

V[3]  =      .298 


/\ 

Data   set   SC 

/ 

> 

a  = 
3  = 

-.303 
.171 

v[$]  = 

V[3]    = 

HI1      - 

.088 
.113 

/     .088 
\-.035 

-.035 
.113 
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X 

3  o 
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Figure  10  3.   Contours  of  log  likelihood 
function  for  data  set  RN 
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Figure  104.   Contours  of  log  likelihood 
function  for  data  set  FL 
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Figure  105.   Contours  of  log  likelihood 
function  for  data  set  SC 


C.   DISCUSSION 

A  recapitulation  of  all  the  parameter  values  is  in 
Table  50.   Some  interesting  observations  to  be  made  from- 
this  table  are: 

(1)  The  slope  for  data  set  RN ,  as  found  by  maximum 
likelihood,  is  much  greater  than  that  of  any  other  method 
or  data  set.   The  first  temptation  is  to  treat  this  as  an 
outlier,  yet  the  evidence  of  the  contour  plot  and  of  the 
validation  of  the  next  section  tend  to  back  up  this  number. 
The  reason  for  this  difference  is  a  possible  subject  for 
further  reserach. 

C2)  Except  for  the  intercept  values,  the  slopes  of  data 
sets  RN  and  FL  seem  to  be  fairly  consistent  within  and  between 
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regression  methods.   This  comment  is  made  in  light  of  the 
difference  of  these  data  sets  and  that  of  SC. 

(3)  Data  sets  RN  and  FL  seem  to  be  similar  in  many  ways, 
yet  data  set  SC  appears  to  be  different  in  both  degree  and 
significance. 

No  other  strong  pattern  is  apparent  in  these  parameter 
values.   Graphical  displays  of  the  parameters,  as  used  with 
grouped  data  are  given  in  Figures  106  (RN) ,  107  (FL) ,  and 
108  (SC) .   Again,  note  the  significant  difference  of  the 
maximum  likelihood  line  for  data  set  RN .   Table  51  contains 
the  data  points  from  which  Figures  106,  107  and  108  were 
drawn. 

After  viewing  these  figures,  the  maximum  likelihood 
approach  is  the  preferred  method  for  the  Peninsula  data  sets, 
whereas  the  robust (C=4)  IRWLS  may  be  best  for  the  Valley  data 
set.   Although  fits  were  made  to  data  set  SC,  it  appears  as 
if  no  great  significance  has  been  found. 


162 


TABLE  50 

PARAMETER  FIT  RECAPITULATION 
FOR  ALL  DATA  SETS 


DATA    SETS 

Method 

1 .      Maximum 

Likelihood 

a 
3 

• 

RN 

.062 
2.918 

FL 

-.171 
.933 

SC 

-  .303 
.171 

2.      OLS 

Group=3 

a 
3 

: 

-    .195 
1.771 

-.146 
1.052 

-  .304 
.506 

Group=4 

a 
3 

: 

-    .170 
1.769 

-.115 
1.357 

-  .228 
.545 

Group=5 

a 
3 

• 

-    .265 
1.719 

-.139 
1.523 

-  .184 
.657 

3.       IRWLS 

Group=3 
C=9 

a 
3 

'• 

.048 
1.013 

-.063 
.664 

-  .116 
.351 

Group=3 
C=4 

a 
3 

j 

.052 
1.031 

-.19  7 
1.049 

.154 
.480 

Group    4 
C=9 

a 

3 

; 

-.103 
.665 

-.093 
.705 

-  .169 
.280 

Group=5 
C=9 

a 
3 

: 

-.065 
.855 

-.107 
.723 

-  .126 
.413 
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MODEL  FITS  TO  GROUP 
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Figure  106.   Estimated  probability  of  greater- than - 
average  total  rest-of-year  rainfall  versus  the 
anomaly  of  logged  rainfall  for  January  for  data  set  RN 
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MODEL  FITS  TO  GROUP=3  FL 
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Figure  10  7.   Estimated  probability  of  greater- 
than-average  total  rest-of-year  rainfall  versus 
the  anomaly  of  logged  rainfall  for  January  for 
data  set  FL 
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MODEL  FITS  TO  GROUP=3  SC 
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Figure  108 .   Estimated  probability  of  greater-than 
average  total  rest-of-year  rainfall  versus  the 
anomaly  of  logged  rainfall  for  January  for  data 
set  FL 
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TABLE  51a 


ACTUAL  VALUES  FOR  MODEL  FITS  OF  RN 


LOGGED 

ACTUAL 

MLE 

IRWLS 

IRWLS 

OLS 

ANOMALY 

VALUE 

C=9 

C=4 

-1.090 

0.000 

.042 

.258 

.255 

.107 

-    .460 

.330 

.218 

.397 

.396 

.267 

-    .190 

.330 

.379 

.464 

.464 

.370 

.010 

.670 

.523 

.515 

.516 

.456 

.16  0 

.670 

.629 

.552 

.554 

.522 

.360 

.670 

.753 

.602 

.604 

.609 

.500 

.670 

.821 

.635 

.638 

.666 

.980 

1.000 

.949 

.739 

.743 

.824 

TABLE  51b 


ACTUAL  VALUES  FOR  MODEL  FITS  OF  FL 


LOGGED 

ACTUAL 

MLE 

IRWLS 

IRWLS 

OLS 

ANOMALY 

VALUE 

C=9 

C=4 

-1.820 

.330 

.12  4 

.219 

.108 

.113 

-1.000 

0.000 

.249 

.326 

.223 

.232 

-    .460 

0.000 

.354 

.409 

.336 

.348 

-    .200 

.670 

.412 

.451 

.400 

.412 

.110 

.330 

.483 

.503 

.480 

.492 

.140 

.330 

.490 

.507 

.487 

.500 

.200 

0.000 

.504 

.517 

.503 

.516 

.340 

.670 

.536 

.541 

.540 

.553 

.500 

.670 

.    .573 

.567 

.581 

.594 

.610 

.670 

.598 

.585 

.609 

.621 

.750 

1.000 

.629 

.607 

.643 

.655 

1.140 

1.000 

.709 

.667 

.731 

.741 
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TABLE  51c 
ACTUAL  VALUES  FOR  MODEL  FITS  OF  SC 


LOGGED 

ACTUAL 

MEL 

IRWLS 

IRWLS 

OLS 

ANOMALY 

VALUE 

C=9 

C=4 

-1.410 

.330 

.367 

.352 

.372 

.266 

-1.070 

0.000 

.381 

.380 

.411 

.300 

-  .960 

.330 

.385 

.389 

.424 

.312 

-  .550 

.330 

.402 

.423 

.473 

.358 

-  .260 

.670 

.414 

.448 

.507 

.393 

-  .100 

.330 

.421 

.462 

.526 

.412 

.020 

.670 

.426 

.473 

.541 

.427 

.200 

0.000 

.433 

.489 

.562 

.449 

.280 

.670 

.437 

.496 

.572 

.460 

.400 

.670 

.442 

.506 

.586 

.475 

.460 

1.000 

.444 

.511 

.593 

.482 

.560 

.670 

.448 

.520 

.604 

.495 

.690 

0.000 

.454 

.532 

.619 

.511 

.780 

.670 

.458 

.539 

.629 

.523 

.970 

.330 

.466 

.556 

.650 

.547 

1.250 

1.000 

.478 

.580 

.680 

.581 
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VII  .    VALIDATION  OF  LOGISTIC  MODELS 

A.   GENERAL 

The  various  parameters  that  were  estimated  in  the  previous 
section  all  may  be  subject  to  some  sort  of  validation.   How- 
ever, this  paper  will  only  view  the  validation  for  the  maximum 
likelihood  approach  on  all  data  sets  and  the  IRWLS  (C=4) 
approach  on  data  set  SC.   The  validation  will  be  conducted 
against  the  reserved,  independent,  data  sets  of  years  19  75 
through  19  80.   These  are  the  same  data  sets  as  used  in  section 
IV. 

Table  52  portrays  the  reserved  data,  in  a  form  for  logistic 
analysis,  and  Figure  109  is  a  display  of  the  derived  contin- 
gency tables  for  the  reserved  data. 
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TABLE  5  2 

RESERVED  DATA  IN  FORM  FOR  THE 
LOGISTIC  ANALYSIS 

DATA  SET  RN 


YEAR 

1975 

1976         -2 

1977 

1978 

1979 

1980 


X 

COMPLEMENT 

Y 

739 

12.92 

0 

446 

11.28 

0 

477 

11.14 

0 

883 

19.62 

1 

542 

18.11 

1 

752 

___ 

— 

DATA  SET  FL 


1975  -.470 

1976  -2.610 

1977  -.470 

1978  .734 

1979  .547 

1980  .633 


DATA  SET  SC 


1975 

-.535 

1976 

-3.773 

1977 

.352 

1978 

1.054 

1979 

.512 

1980 

.405 

12.58 

0 

10.58 

0 

11.10 

0 

19.32 

1 

13.72 

0 

17.97 

1 

11.19 

0 

10.03 

0 

23.38 

1 

16.36 

0 
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Data  set  RN 


Control 


Complement 

+ 
3        0 


+ 

0 

1    2    | 

3 

3 

1   2     | 

5 

Data 

set  FL 

Complement 

- 

+ 

+ 

3 

1   o   1 

3 

Control 

1 

1    1    1 

2 

4    I 

1    1    1 
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Data 

set  SC 

Complement 
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+ 
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1    1 

2 
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+ 
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3 

Figure  109.   2x2  contingency  Tables  of  reserved  data 
controlled  by  the  anomaly  of  January  rainfall.   The 
complement  is  the  anomaly  of  the  rest-of-year  rainfall, 
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B.   RESULTS 

1.   Data  set  RN 

The  model  proposed  by  the  maximum  likelihood 
parameters  is 

. 0618+2. 9183X,. 


,     .  0618+2.  9183X,. 
1  +  e  t 


VIII. 1 


This  model,  when  applied  to  the  reserved  data  yields  Table  5  3 


TABLE  5  3 

RESULTS  OF  LOGSTIC  VALIDATION 
ON  DATA  SET  RN 


YEAR 

X 

Y 

0 

1975 

-.739 

0 

.11 

1976 

-2.746 

0 

.00 

1977 

-.477 

0 

.21 

1978 

.883 

1 

.93 

1979 

.542 

1 

.  84 

1980 

.752 

_ 

.91 

The   9    is  interpreted,  again,  as:   The  conditional 
probability  that  the  complement,  the  total  rainfall  for 
February  through  December,  will  be  above  its  mean  value, 
given  that  the  logged  January  anomaly  was   X ,  .   Thus  it 
appears  that  this  model  tends  to  predict  the  direction  of 
the  complements  deviation  well.   Figure  110  is  a  plot  of  the 
estimated  probabilities  against  the  actual  complement  anomaly 
For  an  acceptable  fit,  this  plot  should  show  an  upward  to  the 
right  slope,  which  it  does. 
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Figure  110.  Plot  of  9   versus  complement 
anomalies  for  data  set  RN 


2.   Data  Set  FL 

This  data  set  is  quite  similar  to  the  RN  data, 
except  that  the  slope  parameter  is  only  a  third  of  that 
of  RN.   The  model  is 


-.171+.9325X. 

e t 

,  .   -.171+.9325X,. 
1  +  e  t 


VII. 2 
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TABLE  54 

RESULTS  OF  LOGISTIC  VALIDATION 
ON  DATA  SET  FL 


Year 

X 

Y 

9 

1975 

-.470 

0 

.35 

1976 

-2.610 

0 

.07 

1977 

-.470 

0 

.35 

1978 

.734 

1 

.63 

1979 

.547 

0 

.58 

1980 

.633 

- 

.60 

A  plot  of  the  probabilities  against  the  complement  anomalies 
is  in  Figure  111. 
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Figure  111.   Plot  of  9,  versus  complement 
anomalies  for  data  set  FL 


This  fit  is  not  as  good  as  that  for  data  set  RN . 
The  outlier,  or  false  prediction  of  19  79  may  not,  however, 
be  far  out  of  line.   The  sparsity  of  points  for  which  the 
complement  anomaly  was  positive  detracts  from  the  validation 
effort. 


3.   Data  Set  SC 


The  maximum  likelihood  model  is 


,,,      -.3025+.  171X,. 
)(!)_  __J~ t 

fc      i  ^   -.3025+. 171X. 
1  +  e  t 


VII.  3 


and  the  IRWLS  model  is 


,„,       .1537+.4799X,. 
,  (2  )     _e t 

t      I  "   .1537+.4799X,. 
1  +  e  t 


VII.  4 


and  the  tabular  results  are  in  Table  55. 


TABLE  55 
RESULTS  OF  VALIDATION  ON  DATA  SET  SC 
Year       X       Y 


1975 
1976 
1977 
1978 
1979 
1980 


-.535 

■3.773 

.352 

1.054 

.512 

.405 


1 
0 
0 
1 
0 


(1) 

9  (2) 

t 

yt 

40 

.47 

28 

.16 

44 

.58 

47 

.66 

45 

.60 

44 

.59 

and  the  plot  of  the  probabilities  versus  the  complement 
anomalies  is  in  Figure  112. 
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Figure  112.   Plot  of  0.  versus  complement 

anomalies  for  maximum  likelihood  and  IRWLS 
parameters  for  data  set  SC 


C.   DISCUSSION 

The  validation  of  the  maximum  likelihood  models  for  data 
sets  RN  and  FL  appears  to  be  acceptable.   However,  data  set 
SC  does  not  appear  to  be  acceptably  modeled.   In  fact,  as 
Figure  112  shows,  the  complement  appears  to  be  almost  inde- 
pendent of  the  control.   This  is  also  shown  by  Figure  105 
where  it  can  be  seen  that  the  contours  are  very  flat  and 
circular  about  the  origin  of  the   (a,S)  coordinate  system. 

The  one  apparent  outlier  of  data  set  FL  may  be  viewed 
as  very  close,  therefore  that  model  can  also  be  assumed 
to  be  validated. 
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VIII.   FURTHER  FINDINGS 

A.   SUMMER  MONTHS 

The  further  investigation  of  the  summer  months,  to 
parallel  the  modeling  of  the  winter  months,  yielded  some 
interesting  results.   These  results  are  shown  here  with 
no  attempt  at  analysis. 

The  summer  months  appear  to  be  increasing  in  total 
rainfall  and  in  variance.   This  is  more  true  for  the 
Peninsula  data  sets  than  for  the  Valley  data  set.   Figures 
113  (RN) ,  114  (FL) ,  and  115  (SC)  show  the  by-month  series 
of  summer  months.   The  total  summer  rainfall  series  by 
year  are  shown  in  Figures  116  (RN) ,  117  (FL) ,  and  118  (SC) . 
The  reserved  data  are  not  included,  yet  it  can  be  shown  to 
continue  the  indicated  trends. 
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Figure  113.   Monthly  plot  of  summer  months 
only,  means  removed,  for  data  set  RN 
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Figure  114.   Monthly  plot  of  summer  months 
only,  means  removed,  for  data  set  FL 
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Figure  115.   Monthly  plot  of  summer  months 
only,  means  removed,  for  data  set  SC 


DATA  IS  FROM  NEWPN 


3.  180 


,060 


/     v  v  \ 


SUMMERS  TOTAL  BY  YEAR 


25 


Figure  116 .   Yearly  plot  of  total  summer 
rainfall  for  data  set  RN 


179 


DATA    IS    FP.vn    NEMFL 


2.310 


<r 

0L 


0.000 


J. 


a  y4v 


\ 


A   VW 


!•■•»» 


SUMMERS    TOTAL    BY     rEAP 


i    '   ■   '   '   I 


40 


Figure    117 


Yearly  plot  of  summer  month 
rainfall  for  data  set  FL 
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Figure  118.   Yearly  plot  of  total  summer 
rainfall  for  data  set  SC 
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A.   SIGNIFICANCE  OF  JANUARY 

The  identification  of  January  as  a  possible  predictor 
for  its  eleven  month  complement  raises  further  questions. 
One  of  the  questions  is  in  determining  which  part  of  the 
complement  lends  the  most  towards  its  predictability. 
Figure  119  is  a  plot  of  the  log-odds  and  chi-square  statis- 
tics for  the  cumulative  complements.   Each  progressive 
column,  to  the  right,  of  Figure  119  indicates  these  statis- 
tics for  another  cumulated  month,  i.e.,  the  first  column 
compares  the  anomalies  of  January  and  February  by  itself, 
the  second  column,  is  a  comparison  of  January  to  February 
plus  March,  and  so  on  until  the  last  column  is  a  comparison 
of  January  to  the  entire  eleven  month  complement. 

Several  occurrences  to  be  noted  from  the  figure  are: 

(1)  The  log-odds  are  consistently  greater  than  zero. 

(2)  The  lack  of  increased  odds  and  significance  during 
the  summer  months. 

(3)  The  similarity  of  RN  to  FL  and  their  combined 
difference  to  SC  in  the  fall. 

These  indications  suggested  a  further  look  at  January 
versus  the  fall  months  only.   This  analysis  is  displayed 
in  Figure  120.   The  vertical  scales  of  Figure  119  and 
Figure  120  are  the  same,  yet  the  horizontal  scales  differ. 
This  figure   has  five  major  divisions.   The  left-most 
division  looks  at  January  versus  singular  months  in  the 
fall.   The  second  division  looks  at  January  versus  pairs 
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of  months  in  the  fall,  and  so  on  until  the  right-most 
column,  which  is  January  versus  the  total  fall  rainfall. 
This  figure  yields  no  apparent  significance,  and  unstable 
odds. 

The  combined  information  of  Figures  119  and  120  are 
mildly  confusing.   One  possible  explanation  may  be  that 
the  summer  months  somehow  cumulate  significance  and  devia- 
tion direction,  in  order  to  allow  the  fall  contribution. 
This  possible  synergistic  affect  should  be  explored  further. 
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Figure  119.   Log-odds  and  significance  versus 
additional  months  cumulated  through  the  year 
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Figure  120.   Log-odds  and  significance 
versus  additional  months  of  the  fall 
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IX.   SUMMARY 

The  analysis  of  rainfall  data  is  carried  out  in  a 
comprehensive  way.   The  autoregressive  Markovian  model 
of  the  early  sections  could  not  stand  up  to  validation, 
but  it  did  point  to  some  sort  of  dichotomy  between  the 
seasons. 

2x2  contingency  analysis  was  effective  in  that  it 
brought  attention  to  the  predictive  ability  of  January. 
This  identification  of  January,  when  followed  by  the 
logistic  analysis  was  seen  to  be  successful  in  two  of 
the  three  data  sets.   Thus,  the  primary  conclusion  of  this 
thesis  is  the  predictive  ability  of  January  rainfall. 

The  physical  reasoning  behind  this  finding  must  be 
left  to  the  meteorologist.   Further  study  of  the  approach 
used  here  may  lead  to  improvement  in  seasonal  or  annual 
rainfall  forecasts  for  certain  climatic  regions. 
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APPENDIX  A 
DATA  SET  RN 


YERR   OCT 
RPR 


NOV 
NAY 


DEC 
JUN 


JflN 

JUL 


FEE 
RUG 


MRR 
SEP 


1951 


.900 
1.  1Q0 


2.  7SG 
.  160 


6.330   10.040 
.240     .060 


2.960 
.  030 


4.410 
.OSO 


1952 


.  170 
1.670 


2.  190 
.490 


5.9O0 
.  220 


2 .  060 
0.  000 


040 
110 


1.200 
.070 


195: 


.330 
.350 


2.  ISO 
.420 


560 
400 


4.260 
0 .  000 


2.  260 
.  160 


4.  910 
.050 


1954 


.030 
1.740 


2.  500 
.370 


3.  130 

.  060 


5.  320 
0.  000 


1  .  740 
0.000 


.  160 
Q .  000 


1955 


.  OSO 
1.670 


1.960 
.590 


9.790 

0.  000 


6.  090 
.  070 


060 
020 


150 
360 


1956 


1 .  000 
1.490 


0.  OOO 
2.390 


340 
200 


4.  650 
0.  OOO 


3.520 
.  070 


1  .920 
.030 


1957 


1.530 
4.710 


.  930 
.560 


3.700 
.350 


3.710 
.040 


5.  660 

0.000 


7.  170 
.480 


1953 


.040 
.290 


.510 
.  120 


.490 
0.  OOO 


4.350 

0  .  OOO 


;.  76o 

.  040 


.320 
3.  140 


1959 


O.  OOO 

.830 


0  .  000 

.340 


.590 
0.000 


4.300 

.030 


4.530 

O.O0O 


.  340 
.  130 


I960 


.070 
1.290 


2.060 
.720 


.350 
0.  0O0 


1 .  390 
0.  OOO 


1.  170 
.  140 


.090 


1961 


.  040 
.300 


1.740 
.  150 


1.  190 
.230 


2.690 
0.  OOO 


i.  170 
.250 


2.570 
.  150 


1962 


1.330 
3.930 


.370 
.  66Q 


2.210 
.040 


3.050 
.040 


2.  70O 
.001 


4.  140 
.  400 


1963 


1.460 
.220 


3.770 

.360 


.530 
.220 


3.500 
.090 


.420 
.350 


2.230 
.010 


1964 


.780 
2.260 


3.290 
.  .  170 


6.450 
.  150 


:.560 
.050 


1  .050 
.  160 


2.  440 
.020 


1965 


.230 
.270 


6.490 
.  130 


5.560 
.  120 


320 
280 


1.380 
.090 


.430 
.320 


196< 


.090 
7.110 


4.740 
.400 


4.  180 

1.560 


5.290 
.020 


450 
060 


5.430 
.  170 


186 


YEAR 

OCT 

NOV 

DEC 

JAN 

FEE 

MAR 

APR 

MAY 

JUN 

JUL 

AUG 

SEP 

1967 

.330 

1.610 

2.270 

3.  100 

1  .  400 

3.  Q60 

.790 

.320 

.  010 

.060 

.  230 

.050 

1963 

.310 

3.  130 

3.  270 

9.450 

7.310 

1.310 

i 
j 

2.700 

.  120 

.420 

.040 

.001 

.  120 

1969 

.500 

.720 

3.030 

5.910 

2.04O 

2.  970 

.350 

.050 

.300 

.030 

.  060 

.020 

1970 

.590 

6.  170 

4.  990 

1  .030 

.620 

1  .  960 

: 

1.  190 

.710 

.030 

.070 

.  130 

.  430 

1971 

.090 

1.990 

4.  760 

1 .  230 

1  .  O50 

.030 

i 

.830 

.090 

.  150 

.060 

.040 

.  1O0 

1972 

2.  460 

5.950 

2.030 

6.050 

5.330 

4.520 

.  130 

.060 

.020 

.  020 

.  050 

.  340 

1973 

2.200 

3.370 

4.730 

3.  730 

.  910 

4.430 

3.400 

.  030 

.370 

.250 

.020 

.010 

1974 

1.540 

.560 

2.480 

1 . 340 , 

3.620 

4.  060 

- 

1.760 

.010 

.  170 

.  170 

.430 

.020 
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Logged  Anomalies 


YEAR   OCT 
APR 


NOV 
MAY 


DEC 
JUN 


JAN 
JUL 


FEB 
AUG 


MAR 
SEP 


1951 


.204 
-.  112 


.237 
-.  1?Q 


.753 
.040 


363 
005 


.  272 
-.017 


.556 
-.  108 


195 


-.231 
.  128 


.  068 

.080 


.620 
.023 


-.421 
-.  054 


1.065 
.010 


-.  344 
-.  1  18 


1953 


-.  116 
-.239 


.055 
.032 


-.367 
.  161 


.  121 
-.054 


073 

054 


.  644 
-.  137 


1954 


-.409 
.  154 


160 
307 


106 
117 


381 
054 


096 
094 


-.934 
-.  185 


195 


•.361 
.  128 


-.007 
.  145 


1.067 
-.  176 


.420 
.014 


1.046 
-.  074 


993 
122 


1956 


.255 
.058 


1.092 
.902 


-.702 
.007 


.  193 

-.054 


.404 
••026 


061 
156 


1957 


.510 
.888 


-.435 
.  126 


236 
125 


.01  1 
-.014 


.792 
-.094 


968 
207 


1953 


-.399 
-.599 


-.680 
-.205 


913 
176 


.223 
-.054 


.307 
-.055 


-.855 
1.235 


1959 


-.433 
-.223 


■1.092 
-.026 


343 
176 


.  129 
.024 


.60S 
-.094 


-.523 
-.  063 


I960 


-.370 
-.025 


026 
224 


697 
176 


-.478 

-.054 


■.  329 
.037 


.  143 
-.099 


1961 


-.399 
-.592 


-.084 
-.  179 


528 
031 


-.233 
-.054 


.715 
.  129 


.  140 

-.046 


1962 


.408 
.741 


-.778 
.  138 


146 
136 


-.  140 
-.014 


.204 
-.093 


.504 
.  151 


1963 


.462 
-.655 


.470 
.302 


387 
023 


-.035 
.032 


-.754 
.206 


.  040 
-.175 


1964 


.  139 
.328 


.364 
-.  162 


696 
036 


-.269 
•.005 


.336 
.054 


.  103 
-.  166 


1965 


■.231 
•.615 


.921 
-.196 


569 
062 


.339 
.  193 


-.046 
-.003 


-.775 
.092 


196c 


-.352 
1.239 


.655 
.018 


333 
764 


.300 
-.034 


- .  733 
-.036 


.736 
-.023 
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YEAR   OCT 
APR 


NOV 
MAY 


DEC 
JUN 


JAN 
JUL 


FEB 
AUG 


MAR 
SEP 


1967 


-.  1  16 
-.272 


-.  133 
-.041 


127 
166 


■•  123 
.805 


•.229 
.  113 


.269 
-.  137 


1963 


-.  163 
.454 


.326 
-.205 


140 
175 


30S 
014 


1.013 
-.093 


-.  295 
-.072 


1969 


-.033 
-.554 


•.550 
•.270 


094 
0S7 


.  394 
-.024 


003 

036 


.246 
.166 


1970 


.026 
-.070 


373 
213 


47S 
146 


-.307 
.014 


622 
023 


-.047 
.  172 


1971 


-.352 
-.223 


.003 
-.232 


439 
036 


737 
005 


336 
055 


•1.  103 
-.090 


1972 


.303 
-.732 


.346 
-.260 


187 
156 


414 

034 


824 
045 


.576 
.  107 


1973 


.725 
.628 


.491 
-.239 


434 
139 


015 
169 


457 
074 


.563 
-.  175 


1974 


.494 
.  161 


-.648 
-.309 


065 
019 


-.639 
.  103 


426 
264 


.439 
-.  166 
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Reserved  Data 


YEAR   OCT  NOV  DEC  JAN  FEE  MAR 

APR  MAY  JUN  JUL  AUG  SEP 

1974  1.54G     .560  2. 430  1.340  3.620  4.060 
1.760     .010  .170     .170     .430  .020 

1975  1.700     .520  .370     .130  2.970  1.520 
1.740     .070  .170     .020     .970  .420 

1976  .600     .720  2.0S0  1.740     .330  1.750 
.049  1.210  .080     .030     .020  .650 

1977  .140     .540  5.350  6.730  4.730  5.240 
5.430     .020  .080     .040     .001  .290 

1973     .020  2.130  1.590  4.820  4.520  4. 410 

.580     .290  .020     .350     .090  .020 

1979    1.800  2.350  3.180  5.950  4.7S0  2.400 

1.770     .570  .040     .730     .090  .090 
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APPENDIX  B 
DATA  SET  FL 


RAW 


YEAR   OCT  NOV  DEC  JAH  FEB  MAR 

APR  MAY  JUN  JUL  AUG  SEP 

193?     .640  1.278  4. 600  4. '590  7.590  5.638 

2.250     .050  0.000     .130  0.000     .170 

1933    1.280     .770  2.090  3.190  1 . 620  2.370 

.280     .630     .070  0.880  0.O00     .138 

1939  .750     .4S8  1.470  3.430  3.208  4.240 
.720      .130  0.000  0.080  0.008  8.000 

1940  0.O08     .790  4.520  6.150  7.730  5. 580 
5.760     .470     .040  0.000  0.000  8.880 

1941  .300  .230  7.870  3.480  2.910  2.438 
3.668  .798  8.888     .120     .188  8.000 

1942  1.130  1.470     .970  3.050  2.750  4.240 
1.300  .100  0.000  0.008     .868  .020 

1943  .620  .420  3.330  3.690  3.088  .710 
1.160  .330     .240     .100     .100  0.088 

1944  1.130  4.658  1.940     .930  3.130  2.330 

.438  .300     .118  8.800     .068  .030 

1945  1.370  1.400  4.230  1.080  2. 560  3.888 

.840  .640  8.O0O  0.8O0  0.00O  .158 

1946  .260  4.310  1.978     .448  1.138  1.958 
.330  .210     .120  0.O00  0.00O  0.000 

194?    1.160  .470  1.940     .098  2.148  4.618 

3.258  .480     .178  0.000  0.000  0.000 

1943    2.130  .350  3.200  1.700  2.340  4.730 

.088  .420     .050     .888     .178  8.888 

1949     .830  1.698  1.370  3.120  1.748  1.640 

1.470  .190     .110  0.00O     .070  .098 

1958    1.768  3.888  2.378  2.368  2.198  1.438 

1.060  .270     .128     .110     .100  .050 

1951  .670  3.350  6.048  3.818  2.118  4.558 
1.188  .128     .178     .840     .070  .190 

1952  .140  2.210  5.040  1.710  0.0O0  .778 
1.258  .578     .228  0.000     .050  .070 
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YEAR 


1953 


1954 


195 


1956 


1957 


1958 


1959 


1960 


1961 


1962 


1963 


1964 


1965 


1966 


1967 


1968 


OCT 

NOV 

DEC 

JAN 

FEB 

MAR 

APR 

MAY 

JUN 

JUL 

AUG 

SEP 

.230 

2.050 

.310 

3.330 

1  .  530 

4.  410 

.780 

.390 

.510 

.  020 

.  160 

.030 

.020 

2.  170 

2.820 

4.  560 

1.560 

.  130 

1.660 

.960 

.050 

.060 

.  100 

.050 

.040 

1  .  650 

8.310 

4.  320 

1.  740 

.  189 

1.560 

.590 

0.  000 

.060 

.  090 

.220 

.650 

0.  009 

.770 

4.  320 

2.610 

1  .S90 

1.340 

2.  640 

.210 

.020 

.  050 

.220 

2.  020 

.850 

3.380 

2.620 

4.  620 

6.760 

3.980 

.490 

.  220 

.  150 

.  100 

.550 

.  040 

.390 

.400 

4.  470 

5.  630 

.  340 

.460 

.  120 

.020 

.030 

.090 

2.  170 

.030 

0.  000 

.620 

3.990 

3.270 

.650 

.370 

.410 

.010 

.030 

.  040 

.010 

.070 

1.450 

.640 

1.340 

1.110 

2.420 

1.200 

.750 

.040 

.040 

.  160 

.  140 

.  100 

2.030 

1.330 

2.230 

5.330 

2.  400 

.270 

.  150 

0.  000 

.  130 

.280 

.080 

1.360 

.440 

2.260 

3.338 

2.270 

4.329 

3.320 

.650 

.050 

.  110 

.  040 

.570 

1.360 

5.260 

.540 

3.190 

.390 

2.30O 

.  100 

.699 

.  190 

.  190 

.330 

.050 

1.060 

3.  160 

5.970 

2.370 

.930 

2.  430 

2.  180 

.  149 

.060 

.060 

.  190 

.010 

.  169 

6.250 

6.230 

3.O50 

1.490 

.510 

.173 

.  140 

.  190 

.320 

.200 

.410 

.  149 

4.690 

3.380 

5.750 

.490 

5.070 

6.319 

.790 

.890 

.060 

.  160 

.230 

.  270 

1.310 

2.960 

3.  110 

1.350 

3.  180 

.720 

.290 

.040 

0.  000 

.2?Q 

.  180 

.430 

3.300 

2.750 

3.430 

7.960 

1.209 

2.490 

.  190 

.260 

.  150 

.240 

.  190 
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> 

YEAR 

OCT 

NOV 

DEC 

JAN 

FEB 

MAR 

APR 

MAY 

JUN 

JUL 

AUG 

SEP 

1969 

.4  10 

1.010 

2.300 

5.  550 

3.  020 

1  .  360 

.430 

.080 

.320 

.050 

.  190 

.040 

1970 

.540 

7.440 

4.450 

.370 

.640 

1.940 

1.  190 

.530 

.050 

.  110 

.  170 

.350 

1971 

.  170 

1.550 

4.  630 

1  .  000 

.  440 

.030 

.720 

.  180 

.290 

.  090 

.  100 

.  100 

1972 

2.540 

5.520 

1.350 

5.550 

6.Q7Q 

3.530 

.310 

.  160 

.030 

.070 

.  220 

.410 

1973 

1 .  960 

4.750 

3.270 

3.  76Q 

1.  190 

4.500 

2.990 

.  110 

.310 

.410 

.21-0 

.  100 
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Logged  Anomalies 


YEAR   OCT  NOV  DEC  JAM  FEE  MAR 

APR  MAY  JUN  JUL  AUG  SEP 

1937     .011  -.137     .454  .332  .332     .703 

.333  -.239  -.130  .095  -.106     .011 

1933     .340  -.436  -.141  .035  -.255     .031 

-.549     .151  -.062  -.071  -.106  -.024 

1939  .876  -.615  -.365  . S46  1.001     .473 
-.254  -.216  -.130  -.071  -.106  -.146 

1940  -.434  -.425     .439  .569  .943     .  63S 
1.115     .047  -.091  -.071  -.106  -.146 

1941  .104  -.760     .914  .102  .145     .049 
.743     .244  -.130  .042  .060  -.146 

1942  .272  -.103  -.591  .001  .103     .473 
.037  -.243  -.130  -.071  -.047  -.126 

1943  -.001  -.656     .306  .143  .988  -.647 
-.026     .266     .085  .024  -.010  -.146 

1944  .272     .725  -.191  -.715  .212     .159 
-.438  -.076  -.025  -.071  -.047  -.117 

1945  .379  -.132     .335  -.666  .051     .203 
-.757     .157  -.130  -.071  -.106  -.006 

1946  -.253     .663  -.131  -1.033  -.462  -.102 
-.511  -.147  -.017  -.071  -.106  -.146 

1947  .286  -.622  -.191  -1.312  -.074     .541 
.651  -.001     .027  -.071  -.106  -.146 

1948  .657  -.727            .166  -.405  .127     .562 
-.719     .013  -.031  .006  .051  -.146 

1949  -.454  -.018  -.406  .013  -.210  -.213 

.108  -.164  -.025  -.071  -.033  -.060 

1950  .532     .379  -.054  -.279  -.058  -.275 
-.073  -.099  -.017  .034  -.010  -.097 

1951  .029     .463     .632  .800  -.084     .530 
-.054  -.225     .027  -.032  -.033     .028 

1952  -.353     .159     .529  -.401  -1.213  -.613 

.015     .113     .069  -.071  -.057  -.073 
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YEAR 


1953 


1954 


1955 


1956 


19? 


1953 


1959 


I960 


1961 


1962 


1963 


1964 


1965 


1966 


1967 


1963 


OCT 

NOV 

DEC 

JAN 

FEB 

MAR 

APR 

MAY 

JUN 

JUL 

AUG 

SEP 

-.237 

.  103 

-.  999 

.063 

-.271 

.505 

-.219 

-.009 

.232 

-.051 

.043 

-.117 

-.464 

.  147 

.071 

.318 

-.278 

-1.061 

.  132 

.335 

-.031 

-.013 

-.010 

-.097 

-.444 

-.032 

.962 

.363 

-.210 

-1.013 

.  144 

.  126 

-.  130 

-.013 

-.019 

.053 

.017 

-1  .007 

- .  698 

.274 

.065 

-.  122 

.054 

.954 

.061 

-.051 

-.057 

.053 

.622 

-.  392 

.208 

-.111 

.503 

.865 

.810 

.061 

.069 

.069 

-.010 

.292 

-.444 

-.673 

-.933 

.301 

.673 

-.391 

-.417 

-.225 

-.  1  10 

-.041 

-.019 

1  .008 

-.454 

-1.007 

-.787 

.210 

.233 

-.633 

-.  170 

.006 

-.  120 

.006 

-.Q66 

-.  136 

-.416 

-.111 

-.774 

-.354 

-.472 

.046 

-.007 

.222 

-.091 

-.032 

.043 

-.  015 

-.388 

.  102 

-.402 

-.225 

.703 

.040 

-.557 

-.  193 

-.  130 

.051 

.  141 

-.069 

.375 

-.642 

-.037 

.  177 

-.034 

.438 

.667 

.  163 

-.081 

.034 

-.066 

.305 

.375 

.827 

-.337 

.035 

-.839 

.010 

-.701 

.  187 

.044 

.  103 

.  180 

-.097 

.239. 

.413 

.673 

-.  183 

-.561 

.049 

.361 

-.207 

-.072 

-.013 

.063 

-.  136 

-.335 

.974 

.709 

.001 

-.306 

-.771 

-.639 

-.207 

.044 

.207 

.077 

.  198 

-.353 

.732 

.316 

.512 

-.820 

.620 

1.260 

.244 

.507 

-.013 

.043 

.061 

-.245 

.026 

-.  151 

.016 

-.364 

.247 

-.254 

-.083 

-.091 

-.071 

.  133 

.019 

-.  126 

.452 

.053 

.851 

.974 

-.395 

.454 

-.  164 

.  101 

.069 

.  110 

.023 
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YEAR   OCT 
APR 


NOV 
MAY 


DEC 
JUN 


JflH 
JUL 


FEE 
AUG 


MAR 
SEP 


1969 


-.140 
-.438 


-.309 
-.261 


066 
148 


.  482 
-.022 


173 
Q6S 


-.  133 

-.  10? 


1970 


-.052 
-.012 


1.  126 
.  120 


42? 
081 


.034 


?24 
051 


-.  105 
.  154 


19?1 


-.32? 
-.254 


-.071 
-.  172 


463 
125 


-.705 
.015 


354 
010 


1.  107 
-.051 


1972 


.730 
-.526 


.868 
-.  189 


222 

100 


.  432 
-.003 


738 
093 


198 


1973 


.602 
.588 


.742 
-.234 


183 
463 


.  162 
.273 


434 
085 


521 
051 


196 


Reserved  Data 


YEAR   OCT  NOV  DEC  JAM  FEE  MAR 

APR  MAY  JUN  JUL  AUG  SEP 

1974  1.9*0     .720  2.430  1.700  2.760  4. 680 
1.530     .110     .180     .330     .450     .080 

1975  1.480     .470     .410     .200  2.820  1.060 
1.660     .070     .130     .140  1.330     .S20 

1976  .110     .730  1.660  1.700     .630  1 . 900 
.030     .960     .110     .020     .040  1.050 

1977  .170     .620  5.570  5.730  4.710  5.160 
5.210     .030     .100     .090     .050     .310 

1973     .100  2.220  1.340  4.700  3. 650  3.710 
.560     .130     .010     .330     .120     .040 

1979    1.580  3.360  3. 830  5.120  4.420  2. 030 
1.520     .580     .040     .900     .140     .130 
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APPENDIX  C 
DATA  SET  SC 


RAW 


YEAR 

OCT 

NOV 

DEC 

JAN 

FEB 

MAR 

FIPR 

MAY 

JUN 

JUL 

MUG 

SEP 

1926 

.420 

S.360 

1.640 

2.  440 

7.530 

1  .390 

1.070 

.  230 

.  100 

0.  000 

0.  000 

.  150 

192? 

.800 

2.  100 

2.330 

.830 

2.430 

3.580 

.960 

.050 

0.  000 

0.  000 

0.000 

0.000 

1923 

.  020 

3.  120 

3.410 

1  .260 

1.370 

2.390 

.  1.290 

0.  000 

1.230 

0.  O00 

0.  000 

0.  000 

1929 

0.  0O0 

0.  000 

.330 

■i  "5761 

3.260 

4.250 

1  .330 

1.610 

0.  000 

0.  000 

0.  ©00 

.030 

1930 

.050 

1.430 

.050 

4.260 

1.560 

.920 

.490 

.  660 

.  400 

0.  000 

0.  000 

0.00  0 

1931 

.020 

1.750 

10.260 

4.320 

4.  690 

.  300 

.  140 

.360 

0.  000 

0.000 

0.000 

0.  000 

1932 

0.000 

.  130 

3.  120 

6.920 

.920 

1.570 

.330 

.320 

.040 

0.  000 

0.  000 

0.  O0O 

1 933 

1.  120 

0.  000 

8.250 

3.  150 

5.  290 

0.  000 

.  130 

.350 

.720 

0.  000 

0.  000 

.060 

1934 

.  150 

2.890 

2.690 

7.  000 

.770 

4.380 

5.910 

0.000 

0.  000 

0.  000 

.650 

0.  000 

1935 

.370 

.810 

2.010 

2.660 

9.970 

1.630 

2.460 

.480 

.340 

.430 

0.  000 

0.000 

1936 

.430 

0.000 

4.  160 

4.620 

6.770 

7.  730 

.520 

0.  000 

.  130 

0.  000 

0.  000 

0.  000 

1937 

.060 

1.010 

6.590 

3.440 

13.020 

3.O90 

2.630 

0.  000 

0.000 

0.  000 

0.  000 

0.  000 

1933 

.310 

1.040 

2.430 

3.540 

2.  700 

3.03O 

.420 

.220 

.  100 

0.  000 

0.  000 

.  100 

1939 

1.040 

.370 

2.380 

9.200 

3.960 

2.320 

.600 

.279 

0.000 

0.  000 

0.  000 

.310 

1940 

1.040 

.370 

2.380 

9.200 

3.960 

2.320 

.600 

.270 

.  0.000 

0.  000 

0.  000 

.310 

1941 

.450 

.338 

9.450 

5.  130 

10.420 

3.920 

5.080 

.610 

.  100 

0.000 

0.  000 

0.  000 
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YEAR 


1942 


1943 


1944 


194; 


1946 


194? 


1943 


1949 


1950 


1951 


1952 


1953 


1954 


1955 


1956 


1957 


OCT 

NOV 

DEC 

JAN 

FEB 

MAR 

APR 

MAY 

JUN 

JUL 

AUG 

SEP 

1.  130 

1.040 

10.49O 

5.710 

2.460 

3.2S0 

4.91Q 

1.010 

0.  000 

0.  000 

0.  000 

0.000 

.740 

1.960 

1.630 

3.  400 

2.  420 

3.  740 

1.070 

0.000 

.090 

0.  000 

0.000 

0.  000 

.740 

.360 

3.  330 

3.940 

9.  230 

1.000 

1.940 

.750 

.320 

0.OOO 

0.000 

0.  OOO 

1.030 

3.240 

2.390 

1.330 

7.610 

4.830 

.280 

.200 

0.  000 

0.  000 

.240 

0.  000 

2.800 

2.050 

7.690 

.850 

2.  770 

3.  280 

0.  000 

.430 

0.  000 

0.  000 

0.  000 

0.000 

.798 

.670 

2.360 

.050 

2.4S0 

4.380 

3.360 

.620 

.  150 

0.  000 

O.  OOO 

O.0O0 

1.740 

.  150 

6.010 

1.  130 

3.040 

5.960 

.  150 

.510 

0.  000 

.020 

.  030 

0.  000 

.  190 

1.470 

1.690 

6.090 

2.970 

2.230 

1.470 

.260 

0.  000 

0.  000 

.050 

0.  000 

2.720 

6.  330 

3.040 

2.950 

1.  940 

1  .040 

1.370 

.690 

.050 

0.  000 

0.000 

.030 

1 .  440 

4.010 

3.360 

9.  800 

1.610 

6.320 

.300 

.220 

0.000 

0.  000 

0.  000 

.070 

.070 

3.000 

3.760 

2.630 

0.  000 

2.090 

1.390 

.620 

.  100 

0.  000 

.090 

0.000 

.300- 

2.  160 

.490 

4.370 

3.280 

4.680 

.530 

.  190 

.250 

0.000 

0.  000 

0.000 

0.000 

1.980 

2.980 

5.530 

1.990 

.310 

2.840 

1.360 

0.  000 

0.  000 

0.  000 

0.  000 

.030 

2.240 

14.940 

6.590 

2.250 

.510 

1.600 

.430 

0.000 

0.000 

0.  000 

.210 

.790 

.020 

.530 

5.210 

5.  1  10 

1.570 

1.720 

1.930 

.090 

8.  000 

0.  000 

.070 

1.470 

.930 

3.990 

4.630 

9.680 

7.320 

8.509 

.560 

.  140 

0.000 

0.000 

.780 
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YEAR   OCT  NOV  DEC  JAN  FEE  MAR 

APR  MAY  JUN  JUL  AUG  SEP 

1953     .030  1.010     .270  6.350  6.920  .230 

.110     .090  O.0O0  O.OO0     .090  2.100 

1959  0.O0O  0.0OO     .710  5.260  6.520  .740 
1.910     .330  0.000  0.O0O  0.OO0  .050 

1960  .050  3.860  1.310  1.910     .910  2.490 
.860     .520     .070  0.000  0.000  0.000 

1961  .170  2.620  1.300  2.000  11.290  2.960 
.110     .340  O.O00  O.O00  0.000  O.OOO 

1962  1.810     .100  2.290  7.810  2.200  4.420 
3.990     .230     .050  O.OOO  O.OOO  .090 

1963  1.30O  3.640     .460  4. 000     .460  3.080 

.390  1.350     .170  0.000     .300  O.OOO 

1364    1.290  3.500  4.330  2. 830  1.050  2. 520 

2.850     .140     .020  O.OOO     .250  O.OOO 

1965  .070  6.630  4.790  1.310  1 . 340  .680 
.270  0.000     .030     .290  O.OOO  .180 

1966  O.OOO  3.740  5.410  6.090     .430  6.160 
6.410     .520     .190  O.OOO  O.OOO  .130 

1967  .270  1.750  1.870  3.360  1.010  3.290 
.670     .430     .120  0.OO0     .020  O.OOO 

1968  .330  1.950  4.O60  15.160  11.970  1.020 
2.140     .120     .100  O.OOO  0.000  0.0O0 

1969  .340     .790  3.110  7.000  4.610  1.510 
.310     .070     .040  0.000  0.000  0.000 

1970  .110  5.360  6.440  1 . 20O     .530  1.490 
1.210     .150  0.000  0.000     .040  .170 

1971  .240  1.430  5.560  1.190     .980  .010 
.890     .030     .020  0.000  0.000  .030 

1972  2.350  5.790  2.660  7.590  9.200  4.430 

.100     .020  0.000  O.OOO  0.000  .070 

1973  1.810  6.310  2.360  4.440  1.930  4.740 
1.970  0.000     .110     .070  0.000  0.000 
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Logged  Anomalies 


YEAR   OCT 
APR 


NOV 
MAY 


DEC 
JUN 


JAN 
JUL 


FEE 
AUG 


MAR 

SEP 


1926 


-.093 
-.090 


1.275 
-.073 


-.437 
.003 


-.347 
-.015 


.705 
.  031 


-.372 
.065 


192; 


.  144 
•.  144 


170 
271 


-.06: 
-.09; 


•.  952 
.015 


-.212 
-.031 


279 
075 


1923 


■.  424 
.011 


455 
320 


.  076 
.710 


-.763 
-.015 


-.531 
-.031 


.  115 
-.075 


1929 


-.444 
.029 


961 
639 


-.304 
-.092 


.269 
-.015 


.005 
-.031 


.415 
-.046 


1930 


-.395 
-.418 


053 
137 


1.359 
.244 


.077 
-.015 


-.504 
-.031 


-.591 
-.075 


1931 


-.424 
-.686 


050 
013 


1.013 
-.092 


.039 
-.015 


.294 
-.031 


-.655 
-.075 


1932 


-.444 
-.532 


796 
279 


003 
053 


.436 
-.015 


— .  792 
-.031 


-.299 
-.075 


1933 


.303 
-.695 


961 
295 


817 
450 


-.  160 
-.015 


.395 
-.031 


-1.243 
-.017 


1934 


-.304 
1.  116 


397 
320 


102 
092 


.496 
-.015 


-.873 
.469 


529 
075 


1935 


-.  129 
.424 


368 
072 


306 
200 


■.235 
.343 


.951 
-.031 


-.257 
-.075 


1936 


-.086 
-.398 


961 
320 


233 
030 


.  143 

-.015 


.606 
-.031 


929 
075 


1937 


-.336 
.472 


263 
320 


619 
092 


-.092 
-.015 


1.  196 
-.031 


.964 
-.075 


1938 


.  150 
-.466 


248 
121 


161 
003 


-.070 
-.015 


-.  136 

-.031 


.  163 
.020 


1939 


.269 
-.347 


646 
031 


190 
092 


.739 
-.015 


.854 
-.031 


-.043 
.  195 


1940 


.269 
-.347 


646 
081 


190 
092 


.739 
-.015 


.854 
-.031 


043 
195 


1941 


.072 
.988 


676 
156 


933 
003 


.238 

••015 


.991 
-.031 


1.052 
-.075 
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YEAR 


1942 


1943 


1944 


1945 


1946 


1947 


1948 


1949 


1950 


1951 


1952 


1953 


1954 


1955 


1956 


1957 


OCT 

NOV 

DEC 

JflN 

FEE 

MAR 

APR 

MAY 

JUN 

JUL 

RUG 

SEP 

.312 

-.248 

1.033 

.321 

-.203 

.211 

.960 

.378 

-.092 

-.015 

-.031 

-.075 

.  1  10 

.  124 

-.441 

.658 

-.215 

.313 

-.090 

-.320 

-.006 

-.015 

-.031 

-.075 

.  110 

-.654 

.057 

.014 

.336 

-.550 

.261 

.239 

.  185 

-.015 

-.031 

-.  075 

.239 

.433 

-.  137 

-.737 

.709 

.529 

-.570 

-.  133 

-.092 

-.015 

.  134 

-.075 

.391 

.  154 

.754 

-.968 

-.  117 

.211 

-.317 

.037 

-.092 

-.015 

-.031 

-.075 

.  13S 

-.443 

-.  196 

-1.534 

-.197 

.440 

.764 

.  162 

.043 

-.015 

-.031 

-.075 

.564 

-.321 

.539 

-.304 

-.043 

.697 

-.677 

.092 

-.092 

.005 

-.002 

-.075 

-.270 

-.057 

-.  419 

.376 

-.066 

-.071 

.087 

-.089 

-.092 

-.015 

.017 

-.075 

.870 

1.031 

-.012 

-.209 

-.  366 

-.530 

.046 

.204 

-.043 

-.015 

-.031 

-.046 

.448 

.650 

.328 

.797 

-.435 

.314 

-.229 

-.  121 

-.092 

-.015 

-.831 

-.007 

-.376 

.425 

.370 

-.280 

-1 . 444 

-.  115 

.244 

.  162 

.003 

-.015 

.055 

-.075 

-.  131 

.  189 

-1.009 

.098 

.010 

.494 

-.392 

-.  146 

.131 

-.015 

-.031 

-.075 

-.  444 

.  131 

-.027 

.293 

-.349 

-.973 

.528 

.538 

-.092 

-.015 

-.031 

-.075 

-.414 

.214 

1.361 

.  444 

-.266 

-.831 

.  138 

.072 

-.092 

-.015 

-.031 

.  116 

.  138 

-.941 

-.933 

.243 

.366 

-.299 

.  183 

.755 

-.006 

-.015 

-.031 

-.007 

.460 

-.304 

.  199 

.  145 

.924 

.934 

1.434 

.  124 

.039 

-.015 

-.031 

.502 
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YEAR 


195S 


1959 


1960 


1961 


196; 


1963 


1964 


196! 


1966 


1967 


1968 


1969 


1970 


1971 


1972 


1973 


OCT   • 

NOV 

DEC 

JAN 

FEB 

MAR 

APR 

MAY 

JUN 

JUL 

AUG 

SEP 

-.414 

-.263 

-1 . 169 

.412 

.625 

-1  .036 

-.713 

-.234 

-.092 

-.015 

.055 

1.056 

-.  444 

-.961 

-.372 

.251 

.573 

-.689 

.251 

-.035 

-.092 

-.015 

-.031 

-.026 

-.395 

.  620 

-.375 

-.515 

-.  797 

.007 

-.  197 

.  098 

-.  025 

-.315 

-.031 

-.075 

-.287 

.325 

-.575 

-.434 

1.064 

.  133 

-.713 

-.023 

-.092 

-.015 

-.031 

-.075 

.589 

-.  366 

-.217 

.593 

-.281 

.447 

.790 

-.  113 

-.043 

-.015 

-.031 

.011 

.3S9 

.574 

-1.030 

.026 

-1 . 066 

.  163 

-.488 

.534 

.065 

-.015 

.231 

-.075 

.385 

.543 

.355 

-.240 

-.726 

.015 

.531 

-.  139 

-.072 

-.015 

.  192 

-.075 

-.376 

1.071 

.343 

-.550 

-.594 

-.724 

-.578 

-.320 

-.015 

.240 

-.031 

.090 

-.444 

.595 

.450 

.376 

-1.037 

.726 

1.  136 

.093 

.082 

-.015 

-.031 

.047 

-.205 

.  050 

-.  354 

-.  1  10 

-.746 

.  213 

-.304 

.037 

.021 

-.015 

-.012 

-.075 

-.  159 

.  121 

.213 

1.200 

1.118 

-.540 

.327 

-.207 

.003 

-.015 

-.031 

-.075 

-.  151 

-.379 

.005 

.496 

.280 

-.323 

-.224 

-.253 

-.053 

-.015 

-.031 

-.075 

-.339 

.965 

.599 

-.794 

-1.019 

-.331 

-.024 

-.  130 

-.092 

-.015 

.003 

.032 

-.  229 

-.073 

.473 

-.799 

-.761 

-1.233 

-.  131 

-.243 

-.072 

-.015 

-.031 

-.046 

.904 

.954 

-.111 

.568 

.373 

.449 

-.722 

-.300 

-.092 

-.015 

-.031 

-.007 

.539 

1.023 

-.  196 

•  111 

-.352 

.504 

.271 

-.320 

.012 

.053 

-.031 

-.075 
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Reserved   Data 


YEAR   OCT  NOV  DEC  JAN  FEB  MAR 

APR  MAY  JUN  JUL  AUG  SEP 

1974  1.340  .760  3.579  2.040  5.430  3.200 
1.660  .020     .040     .070  .120  0.000 

1975  1.310  .730     .330     .030  2.350  1.450 
1.460  .050     .180     .001  1.050     .750 

1976  1.450  .600  1.350  2.450  .460  1.800 

.140  .970     .080  0.000  0.000     .430 

1977  .020  .400  5.630  9.990  3.400  7. 070 
3.600  .030  0.000  0.000  0.000     .360 

1973    O.OOO  2.720  1.200  5.310  4.660  3.920 

.380  .260  0.000     .120  O.OOO  O.OOO 

1979    1.010  2.140  3.370  5.220  10.890  2.390 

2.550  .420  0.000     .950  0.000  0.000 
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